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Abstract 

^Ine  scanning  laser  sensors  which  measure  range  and  reflectance 
are  currently  being  developed  to  be  used  as  the  front  end  of  real-time 
image  processors.  This  paper  establishes  the  coding  combination  of 
intensity  modulation  and  scan  pattern  for  N GaAs  laser  sources  which 
results  in  the  optimum  system  performance.  A noise  model  based  upon 
the  Poisson  point  process  behavior  of  a direct  (optical  intensity)  detector 
is  developed.  Estimation  tlieory.  based  upon  the  Poisson  statistics  of 
the  detector,  is  used  to  develop  the  maximum-likelihood  (ML)  processor 
for  both  range  and  reflectance.  The  realization  of  the  ML  estimator  for 
range  is  shown  to  be  a correlator  or  matched  filter;  a suboptimum  reali- 
zation of  the  reflectance  estimator  is  shown  to  be  a measurement  of  the 
detected  signal  energy,  given  the  range  estimate.  l’erformaix:e  expres- 
sions are  presented  to  relate  the  variance  of  the  range  and  the  reflectance 
to  the  line  scan  system  parameters  and  the  variance  of  both  is  shown  to 
be  inversely  proportional  to  the  detected  signal  energy.  Various  com- 
binations of  spatial  and  temporal  coding  of  the  N laser  sources  to  sepa- 
rate the  return  signal  energy  from  each  ground  resolution  cell  while 
increasing  the  returned  energy  per  cell,  are  evaluated.  Based  upon  tlie 
variaixre  of  the  range  and  reflectance  estimates,  the  optimum  spatial- 
temporal  coding  combinations  is  chosen.  It  is  shown  that  if  a large 
maximum  unambiguous  range  is  not  required,  the  best  performance 
possible  is  achieved  with  N sources  scanned  in  parallel,  each  source 


xiv 


being  sinusoidally  modulated. , If  a large  maximum  unambiguous  range 
is  required,  however,  the  best  performance  is  achieved  with  N sources 
scanned  in  parallel,  each  source  using  sinusoidal  modulation  which  is 
On-Off-Keyed  by  a pseudo-noise  (PN)  code. 
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SPATIAL  AND  TEMPORAL  CODING  OF  GaAs  LASERS 


FOR  A LASER  LINE  SCAN  SENSOR 

I.  Introduction 

The  high  degree  of  mobility  of  modem  tactical  forces  has  generated 
a requirement  for  near -real -time  intelligence.  Traditional  photographic 
techniques  involve  too  many  delays  and  too  much  manual  processing.  Thus, 
there  is  currently  a large  research  and  development  effort  to  develop  a 
real-time  imagery  processing  capability.  A crude  capability  would  at 
least  provide  target  cues  and  coordinate  locations  to  the  photo-interpreters. 
A sophisticated  capability  would  provide  threat  detection  and  target  recogni- 
tion information  to  the  pilot  of  the  collection  vehicle  and  radio  this  same 
information  to  behind-the-lines  intelligence  and  plans  staffs. 

Line  scanning  laser  sensors  are  currently  being  developed  to  be 
used  as  the  front  end  of  real-time  image  processors.  For  some  modes 
of  operation  range  and  reflectance  data  could  be  obtained  from  a laser 
line  scan  sensor  (see  block  diagram,  Fig.  1).  Range  gating  may  result 
in  target  height  profiles  (cross-sectional  boundaries)  which  might  not  be 
readily  apparent  from  reflectance  data  alone  (due  to  a low  contrast  back- 
ground or  camouflage).  The  target  height  and  reflectance  data  could  be 
processed  by  pattern  recognition  algorithms  for  target  cueing  and  threat 
assessment. 
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Fig  1.  b Block  Diagram  of  a Laser  Line  Scan  Sensor 


Fig  1.  Basic  Geometry  and  Block  Diagram  of  a Laser 
Line  Scan  Sensor 
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Future  operational  constraints  require  small,  light,  efficient, 
and  rugged  laser  line  scanners.  A very  favorable  source  candidate  is 


r 


the  semiconductor  GaAs  laser  diode  which  is  very  small  and  light  weight. 

Unfortunately,  the  GaAs  laser  is  also  a low  power  device,  thus  requiring 
the  use  of  many  laser  diodes  to  produce  the  same  performance  achievable 
with  a larger  but  more  powerful  solid  state  or  gas  laser.  A previous 
thesis  (Ref  1)  examined  the  use  of  various  modulation  formats  for 
improved  ranging  and  reflectance  estimation  using  one  GaAs  laser  source. 

This  thesis  will  investigate  various  combinations  of  modulating  and  scan- 
ning techniques  for  N sources. 

Pr(i?lcm 

The  problem  is  to  obtain  the  combination  of  spatial  coding  (scan- 
ning technique)  and  temporal  coding  (signal  modulation)  which  will  result 
in  an  optimum  laser  terrain  mapping  system.  The  purpose  of  the  study 

j 

is  to  arrive  at,  via  parameter  estimation  theory,  a mathematical  model 
for  a theoretically  optimum  three-dimensional  system,  and  then  try  to 
implement  that  model  with  block  diagrams  based  on  available  technology. 

The  performance  of  the  ideal  model  will  serve  as  a standard  or  benchmark 
to  which  the  performance  of  real  systems  may  be  coinpared.  Problems 
encountered  in  implementing  the  tlieoretical  model  will  then  point  to 
areas  requiring  further  research  and  development.  The  goal  is  for  sys- 
tem requirements  to  drive  component  technology,  rather  than  research 
dictating  operational  systems  requirements  and  performance  standards. 
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Scope  and  Assumptions 

The  additive  noise  model  for  the  photodetector  will  be  limited  to 
signal  shot  noise,  dark  current,  thermal  noise  and  background  noise. 

The  background  noise  will  be  assumed  to  be  either  additive  Gaussian 
noise  (for  day  time  performance  estimation)  or  shot  noise  associated 
with  a point  process  (for  night  time  performance  estimation).  Other 
noise  sources,  such  as  random  terrain  reflectance,  random  power 
fluctuations  variations  in  scan  rate,  and  aircraft  vibration  will  be  ignored. 

Initially,  the  sources  will  be  restricted  to  peak-power-limited 
(peak  power  output  can  never  exceed  the  continuous  wave  (CW)  power  out- 
put) GaAs  laser  diodes.  The  impact  of  removing  the  peak-power- limited 
restriction  will  be  addressed  in  a later  chapter.  A recent  study  indicates 
that  some  lasers  may  not  be  peak-power-limited  (Ref  2:  6-14).  GaAs 
laser  diodes  have  too  many  lasing  modes  to  be  practical  in  coherent  sys- 
tems, thus,  only  intensity  modulation  and  direct  detection  (power) 
receivers  will  be  considered. 

The  scanning  process  will  be  modeled  as  if  the  illuminating  laser 
beam  moves  discretely  from  pixel  (ground  resolution  cell)  to  pixel  instead 
of  sliding  along  the  ground  in  a continuous  manner.  It  will  be  assumed 
that  the  dwell  time  (illumination  time)  for  each  pixel  is  so  short  tliat  the 
inter-pixel  smoothing  or  smearing  due  to  tlie  forward  movement  of  the 
scanning  platform  is  negligible.  The  pixel  will  be  assumed  to  be  a 
square  with  sides  equal  to  the  dimensions  of  the  largest  spot  fornted  by 
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die  laser  beam  projected  onto  the  ground  at  the  extremity  of  the  scan. 

Also,  a square  laser  beam  cross-section  will  be  assumed. 

■ Thesis  Overview 

The  first  chapter  of  this  report  will  cover  background  theory. 

The  first  section  will  introduce  a direct  detection  statistical  model  based 
upon  the  statistics  of  the  quantum  nature  of  the  detection  process.  The 
underlying  statistics  of  the  detection  process  will  then  be  used  in  the  i 

following, section  on  parameter  estimation  and  performance  theory.  | 

Here  the  maximum  likelihood  (ML)  estimators  and  Cramer- Rao  lower 
performance  bounds  (variance)  will  be  developed  for  both  range  and  re- 
flectance. The  root- mean- square  (rms)  of  the  Cramer-Rao  bounds 
will  be  the  major  tool  used  in  later  chapters  to  compare  to  performance 
of  various  spatial- temporal  coding  schemes.  This  chapter  will  conclude 
with  comments  on  the  implementation  of  the  ML  estimators  for  range  and 
reflectance. 

The  second  chapter  will  be  devoted  to  the  development  of  potential 
coding- scanning  combinations.  It  will  be  followed  by  a chapter  which 
compares  the  performance  of  the  various  coding- scanning  combinations. 

The  report  will  conclude  with  a summary /interpretation  of  the  analysis 

results  and  recommendations  for  further  research  and  development.  | 
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II.  Theoretical  Background 


The  basic  tools  required  to  compare  the  performance  of  different 
spatial-temporal  coding  combinations  will  be  developed  in  this  chapter. 
The  first  section  will  describe  the  photo-detection  process  and  its  stat- 
istical additive  noise  model.  The  noise  statistics  introduced  in  the  first 
section  will  then  be  used  in  the  following  section  on  parameter  estimation 
and  performance  theory.  The  chapter  will  conclude  with  a discussion  of 
implementation  of  the  estimators  previously  developed. 

Direct  Detector  Statistical  Model 

The  detection  of  optical  energy  is  accomplished  by  converting 
radiant  energy  into  an  electrical  signal  whose  amplitude  is  proportional 
to  the  intensity  of  the  incident  optical  field.  Photons  (dis''rete  energy 
packets,  or  quantum,  of  light)  are  absorbed  by  the  p ho  code  tec  tor.  Some 
of  the  absorbed  photons  excite  conduction  electrons,  thus  creating  a con- 
duction current  i(t)  = q ? 6 (t-t^) , while  other  absorbed  photons  cause 
phonon  (or  acoustic  vibration)  heating.  The  efficiency  (percentage  of 
total  absorbed  photons  which  excited  electrons)  is  given  by  the  quantum 
efficiency  ri.  Since  the  arrival  times  of  photons  and  the  conversion  of 
photons  to  electrons  are  both  probabilistic  in  nature,  the  photodetector 
current  is  a random  process  in  time.  Therefore  the  detector  output 
current  is  a random  process  regardless  of  whether  the  optical  field  is 
stochastic  or  not.  In  addition,  the  output  current  also  has  components 
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due  to  the  irradiance  from  background  objects  within  the  detector  field-of- 
view  and  due  to  thermal  noise  processes  within  the  detector,  all  of  which 
are  stochastic  processes. 

A quantum-mechanically  correct  description  of  the  direct  detection 
process  can  be  given  in  terms  of  the  incident  optical  field  and  the  detector 
generated  current  (Ref  3:  1813“1815,  1831~1834).  Conditioned  on  knowing 
the  signal  component  of  the  complex  envelope  of  the  incident  optical  field, 
the  detector  output  is  conditionally  an  inhomogeneous  filtered  Poisson 
point  process  with  the  rate  (intensity)  function  of 

Ht)  = IUs(5-.t)l’d^  + SjA(,E[|U^(?.t)|Mdf  + Xp  a) 

> / ' ^ 

>-s  (t) 

where 

\ (t)  = rate  function  (electrons,  or  events,  per  second) 

p = Ti/hfjj  = average  number  of  electrons  generated  per  unit 
of  incident  optical  energy 

Ti  = the  quantum  efficiency  of  the  detector 

hfjj  = energy  per  photon 

Ad  = detector  area 

Ug  = complex  envelope  of  the  incident  electromagnetic  field 
due  to  a known  optical  signal 

Ujj  = complex  envelope  of  the  electromagnetic  field  due  to 
background  noise 


Xg(c)  = signal  rate  function 

Xb  = average  rate  function  of  the  optically  filtered  background 
noise 

E = expected  value  (ensemble  average)  operator 
Xp  = dark  current  function 

The  probability  of  detecting  m photons  in  the  time  interval  [0,T]  is  given 
by 

. . m - (i 

PCN(T)  = m]  = ® ; m = 0,l,2,  ...  (2) 

U = X(a)da  (3) 

where  N(T)  is  the  number  of  photons  detected  in  the  interval.  For  Eq  (1) 
to  be  valid,  it  is  required  that  be  such  that  <<  1 where  is  the 

reciprocal  of  the  optical  bandwidth  of  the  background  noise  field.  Condi- 
tioned on  knowing  U„,  the  joint  probability  density  function  of  the  ordered 
event  times  {ti)  and  the  number  of  events  N(T)  in  an  interval  [O,  T]  is 

given  by  rj, 

- So  Ka)da  m 

This  is  a complete  statistical  characterization  of  the  detector  output. 

A partial  description  is  commonly  given  in  the  terms  of  the  first  and 
second  moments  of  the  filtered  detector  output  current.  The  filtering  may 
be  due  to  the  physical  properties  of  the  detector,  or  may  be  any  external 
filter  between  the  detector  and  the  signal  processor.  Let  the  total  detector 
output  current  be 
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The  background  noise  current  may  be  modeled  (Ref  4;  258-259)  as  a 
Poisson  point  process  with  mean  and  covariance  values  of  and 
q*^b  respectively,  so  long  as  1^  (t)  <<  1.  Even  in  the  absence 

of  an  incident  field  upon  the  detector,  there  is  a current  due  to  thermal 
generation.  This  "dark  current"  may  be  modeled  as  a Poisson  point 
process  with  mean  and  covariance  values  of  qX^  and  q^X^  6(t-t').  res- 
pectively. The  last  term  in  Eq  (5),  l'pjj(t).  is  "thermal  noise"  due  to 
losses  in  the  detector  circuit.  It  is  a Gaussian  noise  process  with  zero 
mean  and  correlation  2KTR  6(t-t  0.  where  K is  Boltzmann's  constant, 

T is  the  temperature  in  degrees  kelvin,  and  R^  is  the  "effective"  thermal 
noise  resistance. 

The  additive  noise  model  for  a detector  with  an  output  given  by 
Eq  (6)  is  developed  in  Appendix  A and  is  shown  below  (Fig  2).  The 
noise  process  n(t)  is  zero  mean  with  autocorrelation 
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(6) 


R^a.t  •)=  [q=Xn  + 2^  + q=Xg(t)]  6(t-t') 

Re 

where  ~ ^d'  terms  q^CX^  + Xg(t)]  are  the  autocorrelation 

of  what  Is  commonly  known  as  "shot  noise". 

Thus  n(t)  is  a wiiite,  zero  nxian,  non-stationary,  and  non-Gnussian 
noise  process.  This  means  (strictly  speaking)  that  frequency  dornan  tech- 
niques, i.e. , Fourier  Transform  techniques,  *wliich  require  stationary, 
cannot  be  used  to  calculate  signal -to- noise  ratio  (SNR)  and  other  useful 
performance  parameters.  It  also  means  tliat  all  of  tlie  classical  radar 
results,  such  as  using  SNR  to  calculate  probability  of  error,  are  not  appli- 
cable because  they  are  based  upon  Gaussian  statistics. 

There  are  two  cases  when  frequency  doman  techniques  may  Ixi  used. 
The  first  case  requires  the  IXl  component  of  X„(t),  denoted  X^,  to  be  much 
larger  than  the  time  varying  fluctuations  of  Xg(t).  Then  Fq  (6)  becomes 

Rn^)=  [q"an  + >^o)  + 2^  ^ (2) 

where  t=  t-t'.  Now  n(t)  is  a white,  approxi mately  stationary  Poisson 
noise  process  with  noise  power  spectral  density 

S„(f)  = q=  (X„  + V + 2 (8) 

This  is  commonly  referred  teas  the  "signal  shot  noise  limited  case".  The 
second  case  occurs  when  q°X^  is  tlic  dominant  term  in  Fq  (6).  This  is 
very  often  the  case  for  real  world  detectors  when  tlie  background  noise 
intensity  is  large  and  cooled  detectors  (i  « O)  are  used.  Fqs  (7)  and 
(8)  tlien  become 
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Sn  (0  “ q* 


This  Is  usually  referred  to  as  the  "background  noise  limited  case". 


If  the  rate  function  becomes  very  large,  then  the  filtered  Pbisson 


process  may  be  modeled  as  a Gaussian  process.  In  general  when 


^n  ^o 

liandwicltli  of  h(t)  ‘ 


the  fiker  outfrut  looks'  like  a GaussiSn  random  process  with  mean  IT  = qX(t) 


and  variance  = x(t)  JJh^(t)dt  (Ref  5:  127-131).  Thus,  although  n(t)  is 


not  Gaussian,  since  a filtered  Gaussian  random  process  is  also  Gaussianj- 


n(t)  may  be  modeled  as  being  Gaussian  by  looking  backward  into  the  filter 


from  the  output.  For  this  specific  case,  given  by  Eq  (11),  n(t)  is  a vdute. 


* ♦ <r 

approximately  stationaryT  and  approximately  Gaussian  noise  process  with 


autocorrelation  and  power  spectral  density  given  by  Fqs  (7)  and  (8), 


respectively. 


In  summhry,  the  output  current  of'tlie  photodetector  is  a conditional 


PoissSon  point  process  with  a second  moment  additive  noise  model  as  shown 


in  Fig  1.  The  autocorrelation  of  the«dditive  noise  is  (in  general)  impul- 


sive and  non- stationary.  Only  in  the  less  interesting  cases  of  "signal  shot 


noise  limited",  or  wlicn  the  noise  rate  function  is  mucli  greater  tlian  the 
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bandwidtii  of  h(t),  can  the  noise  be  modeled  as  Gaussian  and  traditioml 


radar  results  be  used  to  analyze  die  optical  system.  In  the  more  prac- 


tical case  of  "background  noise  limited"  new  performance  techniques  must 


be  dcvelqied  using  Poisson  statistics. 
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r 

{^ranxiter  Rstiriiati on  and  rbrfortnance  Theor 

This  section  will  have  three  parts.  Tlie  first  will  develop  the 
maximum  likelihood  (ML)  estimators  for  both  range  and  reflectance  using 
Paisson  statistics.  It  will  also  discuss  tlie  structure  imposed  upon  the 
signal  processor  by  the  ML  estimators.  The  second  part  will  use  Poisson 
statistics  to  derive  the  Cramer-Rao  (CR)  performance  lower  bounds  in 
general,  and  will  look  at  the  CR  bounds  for  several  specific  modulation 
techniques.  The  third  part  will  compare  the  CR  bounds  developed  using 
Poisson  statistics  with  those  developed  previously  using  Gaussian 
statistics  (Ref  1:  21-56). 

Maximum  Likelihood  Parameter  Estimation  . The  principle 
behind  the  measurement  of  range  and  reflectance  by  the  laser  line  scan 
system  is  very  similar  to  that  of  a classical  radar.  As  with  a radar,  the 
time  delay  (or  phase  difference)  between  a transmitted  signal  and  its  return 
is  proportional  to  the  target  range.  In  addition,  the  return  signal  ampli- 
tude will  be  proportional  to  the  target  reflectance. 

All  measurements  of  the  actual  values  of  the  delay  rand  the  re- 
flectance p will  be  degraded  by  system  noise,  and  thus  will  only  be  esti- 

t>k 

mates  of  the  true  values.  These  estimates,  denoted  f and  P , are  based 
upon  the  observation  of  die  ordered  event  times  {t^)  of  detected  electrons 
and  the  number  of  events  N(T)  in  an  interval  CO,  Tl.  The  desire  is  to 
use  statistical  tools  available  to  determine  a signal  processor  such  that 
E[t]=  T (zero  bias)  and  var  Ct-^]  (a  measure  of  how  close  the  estimate  ^ 
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Is  to  the  real  value  t ) is  small  (witli  similar  expressions  for  P and  p). 
Since  no  a priori  information  about  the  random  parameters  t and  p is 
known,  maximum  likelihood  (ML)  estimation  will  be  used.  Given  die  ob- 
servation A,  the  ML  estimate  of  the  random  parameter  A is  defined  by 
the  likelihood  function 
Maximize 

by  Choice  ^ |A  I ^ M ^ (^2) 

of  A I A = 

Ref  6:  63),  or  equivalently,  the  log  likelihood  function 
Maximize 

by  Choice  { ^ f-  i I A)  j } (13) 

of  A ZIA  = 

where  f^  | ^(Z  | A)  is  the  conditional  probability  density  function  of  (pdf) 
of  the  observation  Z,  conditioned  on  knowing  the  parameter  A.  For  con- 
venience, the  subscripts  will  be  dropped  from  the  traditional  notation  for 
pdf's  and  Nj  will  be  used  to  represent  N(T)  = m.  From  Eq  (13)  it  follows 
tliat  a sufficient,  though  not  necessary,  condition  for  ML  estimate  can  be 
found  from  the  log  likelihood  equation  I (A)  given  by 
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The  rate  function  in  Eq  (14)  is  given  by  rewriting  Eq  (1)  as 


X.(t)  = Mt)  + >4l  (16) 

^s(0  = Pst  = ppP^mCt)  ; |m(t)|£l  (17) 

where  m(t)  is  the  intensity  modulation  impressed  upon  the  transmitter, 
and  the  detector  power  is  the  peak  power  limited  (CW)  laser  output 
Pj  times  all  other  gain  and  loss  mechanisms  such  as  atmospheric  absorp- 
tion, transmitter  and  receiver  field-of-view,  and  receiver  aperture  area. 
Substituting  Eqs  (15)  and  (16)  into-Eq  (13)  yields 

m 

Cn[f  {ti  } , N-p  I T,  p)]  = 07i[ps(t|-T)  + ] (18) 

T 

" Jo  Cps(a-T)  + X^]  da 


The  simultaneous  estimators  t and  p may  be  found  by  substituting  Eq  (18) 
into  the  simultaneous  equations 


a0n[f(  {tj}  , Nx  I T,  p)  ] 

3 T 

a0n[f(  {q  } , Nx  I T,  p)  ] 


= 0 


= 0 


and  solving  simultaneously  for  rand  p.  Using 


(19) 

(20) 


TV  ® ^ ® ^ ) <21) 

substituting  Eq  (18)  into  Eq  (19),  rearranging  and  canceling  like  terms 
gives 
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(22) 


? .S(t-T) 

i=l  Ps(t-T)  + 


3 s(a-T)da 

- Jo  -sr 


The  righthand  side  of  Eq  (22)  becomes  s(T-t)  - s(-'t)  = 0 for  periodic 
modulation  signals,  or  if  the  observation  interval  [0,T]  is  much  greater 
than  delay  t so  that  "end  effects"  may  be  ignored, 


m 


V s(t-T) 

R ps(t-T)  + 


0 


(23) 


t=^ 
T=  T 


Substituting  Eq  (18)  into  Eq  (20)  and  rearranging  terms  gives 


m 

E 


s(t-T) 


i=l  p^s(t-T)  + 


t = ti 

A 

P = P 


= Jo  s(a-T)  aa  = Q,  a constant  (24) 


The  simultaneous  solution  of  Eqs  (23)  and  (24)  for  ^ and  ^ would 
give  the  joint  estimates  for  delay  and  reflectance  respectively.  Unfort- 
unately, Eqs  (23)  and  (24)  are  inseparable  and  there  is  no  solution  for  f 
Independent  of  ? and  vice-versa. 

However,  the  form  of  Eq  (18)  is  that  of  a correlator.  This  may 
be  seen  by  recalling  that  i(t)  = q^?  6(t-tj).  If  i(t)  is  multiplied  by 
q“'cn  X(t)  and  integrated  using  the  sifting  property  of  the  impulse  function, 
the  result  is 


m 


-T 

Jo  ®(t-ti)  q ' en  X(t)  dt  = Ih  X(ti) 


(25) 


Now  subtracting  a constant  Q = Jj’x(a)da  and  substuting  Eq  (16),  the 
result  is  Eq  (18) 
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m 


Cf  ( { ti  3.  Nt,|t,  p)  ] = Z en  [ Ps(t.-T)  + X ] 

^ i=  I ' " 

- Jo  Cps(a-T)  + Xplda 


(18) 


Fig.  3 mxn  Bank  of  Correlators  Implementation  of 
ML  Estimators 


which  may  be  realized  approximately  by  the  mxn  bank  correlations  shown 
above  (Fig  3).  The  "biggest  picker"  picks  the  output  of  the  correlator 
with  the  largest  amplitude  and  the  pj  and  Tj  that  corresponds  to  that 
correlator  are  tlie  estimates  p and  respectively.  The  correlator 
realization  shown  above  requires  mxn  correlators  and  is  not  practical  to 
Implement. 
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The  correlator  input  q"'  ^ [ps(t)  + X ^^3  may  be  written  as 


If  Xjj  >>  ps(t),  then  the  second  term  on  the  righthand  side  of  Eq  (26)  may 
be  expanded  using  a series  expansion.  Keeping  only  die  first  term  of  the 
expansion,  Eq  (26)  equals 


(27) 


The  likelihood  equation,  given  by  Eq  (18),  may  now  be  rewritten  as 


Maximize  [ f J i(t)  1.  p ^ dt  - rI(ps(a-T)  + X„)da]  (28) 

on  P.  T o q ^ 


Keeping  only  the  p and  t dependent  terms,  Eq  (26)  becomes 

Maximize  [ P s(£)  dt  - To  Ps(a-T)  da]  (29) 

on  p,  T ‘'o  q Xn 

The  first  term  in  Eq  (29)  describes  the  correlator.  Now  the  detector 
current  i(t)  is  correlated  with  a linear  function  of  s(t),  as  in  familiar 
radar  applications.  If  the  observation  interval  [O,  T]  is  much  greater 
than  the  largest  expected  delay  t , i.  e. , if  end  effects  can  be  ignored, 
then  the  last  term  in  the  brackets  can  be  ignored. 

A practical  (suboptimum)  estimator  realization  is  shown  below 
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Fig  4.  Suboptimum  ML  Estimators 

The  delay  between  the  transmission  of  the  signal  and  the  occurrence  of 
a peak  is  equal  to  The  amplitude  of  the  peak  detector  output,  given 
the  delay  estimate  is  proportional  to  the  reflectance  estimate  p. 

Note  that,  in  general,  the  detector  output  is  not  correlated  with 
s(t),  as  would  be  done  in  a conventional  radar,  but  instead  is  correlated 
with  the  weighted  natural  log  of  ps(t)  + Thus,  once  again,  it  must 
be  stressed  that  classical  radar  techniques  do  not  always  apply  to  optical 
systems  involving  direct  (power)  detection.  Only  in  the  background 
noise  limited  case  may  the  familiar  radar  ML  estimator  be  used  opti- 
m umly. 

Cramer-Rao  Lower  Performance  Bounds  . It  was  stated  in  the 
previous  section  that  it  is  desirable  that  var  [ t-t]  and  var  [p-p]  be 
small.  These  variances  are  very  difficult  to  compute;  therefore,  a 
lower  bound  for  the  variances  will  be  computed  using  the  Cramer-Rao 
(CR)  lower  bound.  The  CR  lower  bound  for  the  variances  of  the  esti- 
mator errors  for  t and  p is  defined  as 


!1 


y 

£ 

> 
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(30) 


j 


i 


1 


oj  = var  [aj([tj3,  N^)  - A^] 


where  Jjj  is  the  jjth  element  of  the  2x2  square  matrix  (Ref  6:  79). 
The  elements  of  J are 


^ 3*^f({ti},  N^tIt.p 

" BAjSAk 


(31) 


where  A^  = t and  Ag  = P.  The  CR  lower  bound  is  an  optimistic  perform- 
ance measure.  It  indicates  tliat,  regardless  of  how  good  tlie  signal  pro- 
cessor, the  smallest  error  variance  possible  occurs  when  Eq  (30)  achieves 
equality.  Such  a processor  is  called  an  efficient  estimator.  Any  sub- 
optimum process  will  result  in  an  error  variance  larger  than  die  lowest 
limit  given  by  Eq  (30). 

The  elements  of  ^ may  be  evaluated  by  substituting  Eq  (18)  into 
Eq  (31),  using  Eq  (21)  and  die  observation  diat 


E Ps(ti-T)]  = E C ps(t-T)dt]  = Ps(t-T)  dt  (32) 

m 

where  i(t)  = 9 ® E [i(t)l  = (Ps(t-T)  + X^)q  =1(0.  The  elements 

of  J[  are  then 


li 

P=S=(t-T) 

dt 

(33) 

ps(t-T)  +T„ 

• 

II 

S(t-T)s(t-T) 
Ps(t-T)  + 

dt 

(34) 

J„=  So 

S=(t-T) 

dt 

(35) 

Ps(t-T)  + \^ 
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It  can  be  shown  that,  in  tlie  limiting  cases  of  \^  >>  Ps(t-T)  or  « Ps(t-T) 
Eq  (34)  is  equal  to  zero  so  long  as  s(t)  is  periodic  or  the  observation 
interval  T » t.  Therefore,  substituting  Eqs  (33)  and  (35)  into  Eq  (30)  gives 

'T' 


,T  s^(t-T) 


(Background  Noise  Limited  Case); 

‘'t  p^E?"  “ % 

Sbpf  s 

o>  a Jjl 
P Es 

(ii)  for  X^  <<  Ps(t)  (Signal  Shot  Noise  Limited  Case): 
k 


dt]-' 

(36) 

dt]“" 

(37) 

two  equations  become  (i)  for 

X>>  ps(t) 

• 

(38) 

(39) 

r*  2 


P^  E=  b; 

^ SBPF  ® 


jTs(t)dt 


(40) 

(41) 


(Ref  5:  308-309),  where 

^ jr«  u)’'  S(u))  du) 


B3  M 

s " 


= 4Tr® 


j;,  s(i»)d«) 


jr.f° 

r.  iFsWi^ 


E ^ jJs=(t)dt=j:.  lF3(f)l^df 


(42) 


(43) 
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and  F„(f)  is  the  Fourier  Transform  of  s(t).  The  subscript  BPF  denotes 
"bandpass  filtered",  k is  a constant  of  proportionality,  arbitrarily  set 
equal  to  one,  and  has  units  of  seconds  ' . The  parameter  Bg  is  the  nor- 
malized effective  (or  rms)  bandwidth  of  s(t)  where  s(t)  is  defined  by  Eq 
(17).  The  parameter  E is  proportional  to  the  energy  of  the  envelope  of 
the  intensity  of  the  incident  optical  field.  One  may  also  think  of  Eg  as  the 
energy  of  ig(t)  s=  qXg(t),  i.e. , the  post  detector  signal  energy.  Also  note 
that  s(t)dt  is  just  the  average  number  of  photons  detected  over  an  obser- 
vation interval  [0,T]  and  is  proportional  to  the  energy  of  the  complex  enve- 
lope of  the  incident  optical  field.  Thus,  in  both  limiting  cases,  the  mean 
square  error  of  both  f and  $ may  be  reduced  by  increasing  Eg,  and  con- 
sequently jjs(t)dt.  Also,  the  mean  square  error  of  “V  depends  upon  the 
mean  square  bandwidth  of  s(t).  The  mean  square  error  of  t will  be  mini- 
mized by  making  the  mean  square  bandwidtli  as  large  as  possible. 

There  is  another  interpretation  of  Eqs  (36)  and  (37)  for  the  signal 
shot  noise  limited  case  which  will  be  presented  for  the  sake  of  complete- 
ness (Ref  7:  85-86).  If  the  intensity  modulated  signal  s(t)  is  represented 
In  terms  of  the  complex  field  Ug(r,  t)» 

Xg(t)  = PPAj|Ug(?,t)I=  (44) 


For  Xn<<s(t),  Eqs  (35)  and  (36)  become 


1 

4ppAdEiulB‘ju| 


2 


(45) 


(46) 


r»  a 


p PAdEpi 


E |y  j and  B*  j are  the  energy  and  the  mean  square  bandwidth  of  the  inci- 
dent optical  field,  respec  ’ /ely. 

In  the  following  section  (Additional  Comments  on  Implementation) 
it  will  be  shown  that  the  mean  square  bandwidth  and  energy  of  all  pulse  code 
modulation  schemes  to  be  considered  in  this  report  are  approximately 


^PC  ^ ^ » P ^ ^ 1 

P^PcfTd  _ P^PcfPto 
Eg  = — 4 


(47) 

(48) 


where  tQ  is  the  pulse  duration  = pt^  is  the  dwell  time  and  code  word 
length,  and  p is  the  number  of  chips  in  the  pulse  code  sequence  or  word. 
Since  the  slant  range  to  a pixel  is  given  by  r = cr/2  where  c is  the  speed 
of  light  and  t is  tlie  round  trip  time  delay,  it  can  be  shown  tliat  the  range 
error  variance  is  related  to  the  time  error  variance  by 


(49) 


Substituting  Bqs  (47)  - (49)  into  Eqs  (38)  - (41)  gives  the  CR  lower  per- 
formance bound  for  pulse  code  modulation: 
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(53) 


Mean  square  bandwidth  has  little  meaning  for  sinusoidal  (harmonic) 


modulation.  Therefore,  the  CR  bounds  must  be  evaluated  by  letting  s(t) 
be  given  by 


Pd 

s(t)  = p-j"  [ 1 + sin  2tt  f j^t  ] 


(54) 


and  substituting  directly  into  Eqs  (36)  - (37).  After  some  calculus  and 
using  Eq  (49),  the  CR  lower  performance  bounds  for  sinusoidal  modulation 


are: 
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The  CR  bounds  expressed  by  Eqs  (5Q  - (53)  and  Eqs  (55)  - (58)  will  be 
used  extensively  in  Chapter  III  comparing  the  various  spatial-temporal 
coding  combinations. 


At  this  point  it  should  be  noted  that  if  the  peak  power  limited 


restriction  is  removed,  and  if  Pj.  increases  by  the  same  factor  for  both 

PCM  and  sinusoidal  modulation,  dien  E will  increase  the  same  way  for 

both  modulation  types.  Furthermore,  if  therms  power  transmitted  is 

allowed  to  be  a constant  regardless  of  the  modulation  used,  then  E_ 

s 

will  be  identical  for  both  modulation  types.  Eqs  (35)  - (37)  indicate 
that  if  Eg  is  always  the  same,  then  the  reflectance  performance  will  be 
independent  of  the  type  modulation  used,  and  the  range  performance  will 
depend  only  upon  the  mean  square  bandwidth  of  the  modulation  scheme. 

CR  Bounds;  Gaussian  Versus  Poisson  Statistics.  One  might  ask 


how  the  Poisson  performance  bounds  compare  with  those  developed  by 
Chapuran  using  Gaussian  statistics  (Ref  1:  25-26,  54).  In  order  to  make 
a comparison,  the  Gaussian  noise  power  spectral  density  is  defined  to  be 


No  = 


Xq  ; Signal  Shot  Noise  Limited  Case 
; Background  Noise  Limited  Case 


wliere  Xq  is  the  direct  current  (DC)  component  of  Xg(t)  as  defined  by 
r£q  (17).  It  must  be  noted  that  tlie  definition  of  the  mean  square  band- 
width for  pulse  code  modulation  used  by  Chapuran  is  2tt  larger  than  that 
used  in  this  paper  and  that  reflectance  is  hidden  in  Chapuran's  definition 


of  detected  power.  After  making  the  proper  definition  adjustiT)ents,  a 
direct  comparison  of  the  Gaussian  and  Poisson  results  for  both  pulse  code 
modulation  and  sinusoidal  modulation  reveals  that  in  every  case  but  one 
the  rms  error  is  larger  for  the  Gaussian  case.  The  one  exception  occurs 
for  the  rms  error  of  range  estimation  for  the  background  noise  limited 
case  using  pulse  code  modulation.  For  this  specific  case  the  rms  errors 
were  the  same  for  both  the  Gaussian  and  tlie  Poisson  results.  Thus, 
"quick  and  dirty"  performance  "questi mates"  may  be  made  using  the 
appropriate  No/2  for  the  Gaussian  statistical  CR  bounds  so  long  as  it  is 
realized  the  performance  results  will  be  pessimistic.  However,  for  a 
more  exact  and  careful  performance  analysis,  the  CR  performance 
bounds  developed  using  Poisson  statistics  must  be  used  wlienever  a direct 
detection  optical  detector  is  used  in  the  background  noise  limited  regime. 

Additional  Comments  on  Implementation 

This  section  vslll  discuss  rlie  implementation  of  ML  estimators  for 
T and  p utilizing  matched  filters.  Matched  filters  and  dieir  performance 
characteristics  (in  general),  one-shot  matched  filters,  and  closed  loop 
estimators  will  be  discussed  in  turn.  The  section  will  conclude  with  a 
few  brief  comments  on  transmitter  and  receiver  optics,  background  noise, 
and  typical  limiting  cases  for  o^xirational  deployment. 
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1 

f. 

I Matched  Filters  in  General  . The  ML  estimator  for  a direct 

[ 

I 

detector  (developed  earlier  in  this  chapter)  is  a correlator.  It  can  be 
shown  that  correlators  and  matched  filters  are  equivalent  methods  for 
synthesizing  the  same  ML  estimator  (Ref  8:  315-317).  If  the  output  of 
a receiver  y(t)  = m(t)  + n(t)  is  passed  through  a matched  filter  with 
impulse  response  h(t)  = s(tQ  -t),  the  output  of  the  filter  is  just  the  time 
autocorrelation  of  the  signal,  and  is  given  by 

R(''’)  = s(a)s(a+T)da  (61) 

■ ^o 

The  output  of  the  matched  filter  peaks  when  t=  tQ  and  its  amplitude  is 
proportional  to  the  received  signal  energy.  Thus,  the  output  of  a direct 
detector  may  be  processed  via  a matched  filter.  The  output  may  be 
routed  to  a peak  detector.  The  elapsed  time  between  the  transmission 

of  the  signal  and  the  occurrence  of  a peak  is  the  delay  estimate,  and  the 

j 

amplitude  of  the  peak  is  proportional  to  the  reflectance  estimate. 

I The  matched  filter  realization  of  the  MLestiimtor  is  shown 

below  (Fig  5).  Note  that  the  filter  impulse  response  is  the  time  inverse 
' of  the  correlator  reference  input.  Therefore,  the  restriction  that 

applied  to  using  established  radar  techniques  for  optical  detection  using 
correlators  also  hold  for  optical  detection  using  matched  filters. 
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Fig  5.  Matched  Filter  Implementation  of  ML  Estimators 

In  order  to  make  the  signal  s(t)  and  its  autocorrelation  straight 
forward,  several  simplifying  assumptions  will  be  made.  It  will  be 
assumed  that  both  range  and  reflectance  are  constant  in  a pixel.  Also, 
it  will  be  assumed  that  there  is  no  significant  spatial  filtering  caused  by 
scanning  a beam  over  the  ground.  Recall  that  die  scanning  process  is  to 
be  modeled  as  if  the  illuminating  beam  moves  discretely  from  pixel  to 
pixel. 

The  idealized  autocorrelation  of  a pulse  is  a triangular  shaped 
lobe  (Fig  6).  The  CR  lower  bound  assumes  that  the  matched  filter  out- 
put is  a single  lobe,  and  it  computes  the  error  bound  of  noise  perturbations 
about  this  lobe.  If  the  matched  filter  output  (autocorrelation)  has  large 


H to  h- 


Matched 

Filter 


Fig  6.  Matclied  Filter  Response  for  a Single  Pulse 
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sidelobes,  then  a noise  spike  added  to  a sidelobe  could  be  larger  than  the 
peak  of  the  mainlobe,  thus  causing  an  additional  error  in  the  delay  measure- 
ment (Fig  7).  This  is  the  problem  of  global  accuracy  versus  local  accu- 
racy which  is  discussed  in  the  literature  (Ref  9;  294-307).  Therefore, 


Fig  7.  Error  Due  to  Large  Sideldaes  Plus  a Noise  Spike 


it  is  desirable  to  use  pulse  codes  for  s(t)  that  have  bi-level  autocorrela- 
tions, i.e. , whose  autocorrelation  peaks  when  t = 0,  and  is  a constant 
for  T/0  (Fig  6).  Pseudo-noise  (PN)  codes.  Barker  codes.  Totally 
Orthogonal  Complementary  Pair  codes,  and  others  have  bi -level  auto- 
correlations. 

As  stated  earlier,  tlie  range  is  derived  from  the  delay  estimate 
through  the  expression  r = ct/2.  However,  since  the  transmitted  sig- 
.'lal  is  periodic  with  period  T,  a signal  return  for  s(t)  looks  just  like  a 
return  for  s(t  + nT)  where  n is  an  integer.  This  implies 


c(t  modulo  T) 
2 
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Thus,  the  maximum  unambiguous  range  that  may  be  measured  is  given  by 

CT  (63) 

2 

where  T is  the  dwell  time.  From  Eqs  (50)  and  (52),  the  range  perform- 
ance will  be  best  when  T = Tq  = pt^,  where  tQ  is  as  small  as  permitted  by 
the  limiting  bandwidth  of  the  transmitter  or  detector,  whichever  is  smallest. 

Range  resolution  is  a measure  of  how  small  the  slant  range  between 
two  objects  may  be  and  still  be  distinguished  from  one  another  by  the  range 
estimator.  Since  each  pixel  is  assumed  to  have  a constant  range,  range 
resolution  is  of  no  significance  to  this  system  and  will  not  be  evaluated  for 
each  different  modulation  scheme. 


One-Shot  IVfatched  Filters.  The  autocorrelation  of  a signal  as 
defined  by  Eq  (61)  may  be  accomplished  by  "sliding"  the  signal  past  a 
replica  of  itself  while  at  the  same  time  multiplying  and  integrating.  This 
is  usually  accomplished  by  "sliding"  the  signal  through  a matched  filter, 
normally  a delay  line  summer.  Since  the  performance  is  best  when  a 
code  word  is  long  (as  shown  later),  a code  word  reflected  by  a given 
pixel  gets  to  "slide"  through  the  matched  filter  only  once;  hence,  the 
name  "one-shot"  matched  filter.  This  is  the  type  of  matched  filter  that 
will  be  used  to  process  the  pulse  code  modulation  schemes  considered  in 
this  report. 
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PN  codes  are  well  known  and  are  easy  to  generate  using  shift 
registers  (Ref  10;  143-152).  The  normalized  autocorrelation  of  a PN 
code  p chips  long  is  bi-level,  periodic  (or  circular),  and  has  minimum 
correlation  of  -1/p  (Fig  8).  However,  in  order  for  the  autocorrelation 
to  be  bi-level  with  no  sidelobes,  many  PN  codewords  must  slide  tlirough 
a matched  filter  which  is  matched  to  several  codewords.  If  two  adjacent 
pixels  are  more  than  cto  distance  apart  the  string  of  continuous  code- 
words will  be  broken  and  sidelobes  will  appear  in  the  matched  filter 
output.  Thus,  the  chip  duration  t^  should  be  longer  than  the  maximum 
AT  expected  between  adjacent  pixels.  Unfortunately,  as  shown  by  Eq  (50), 
this  degrades  the  range  estimation  performance.  Another  approach 
might  be  to  transmit  at  least  three  code  words  per  pixel.  This  would 
require  that  p'=p/3,  and  again  Eq  (50)  indicates  that  the  performance 
will  be  degraded.  The  choice  of  which  of  the  two  approaches  is  best 
depends  upon  the  maximum  At  expected  and  the  minimum  t^  the  system 
is  capable  of  achieving.  It  should  be  pointed  out  there  are  other  classes 
of  periodic  pulse  codes  that,  although  tliey  do  not  have  a bi-level  auto- 
correlation function,  have  relatively  small  sidelobe  structures  and  may 
be  used  for  ranging.  Many  times  these  codes  are  found  via  exhaustive 
computer  searches,  and  thus  are  not  considered  here. 
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1 


Fig  8.  Normalized  Autocorrelation  of  PN  Code  p Chips  Long 

The  approximate  mean  square  bandwidth  of  a PN  code  was  calculated  by 
Chapuran  (Ref  1:  83-87)  using  the  definition  for  given  by  Eq  (42). 

The  power  spectral  density  S(uj)  for  a PN  code  has  a (sin  x/x)®  distri- 
bution, and,  strictly  speaking,  has  infinite  bandwidth.  Chapuran 
assumed  that  the  significant  portion  of  the  spectrum  envelope  only  ex- 
tended out  to  the  second  zero  crossing  on  each  side  of  the  double  sided 
spectrum.  Thus,  the  resulting  mean  square  bandwidth  is  only  approxi- 
mate. However,  it  is  good  for  most  realizable  systems  since  real 
systems  do  not  generate  pulses  with  zero  rise  and  fall  times.  Thus 
the  results  are  approximate  for  a wide  variety  of  pulse  formats,  not 
just  for  square  pulses.  Unfortunately,  it  can  be  shown  that  the  definition 


which  he  used  for  the  power  spectral  density  for  a PN  code  Spj^(iu)  was 
2tt  times  too  small.  Making  tlie  proper  adjustment,  the  mean  square 
bandwidth  for  a PN  code  is  given  by 


(64) 


1 


The  output  power  of  the  laser  can  be  represented  by 


peak  with  no  sidelobes.  Therefore,  CP  codes  do  not  exhibit  the  problem 
encountered  with  PN  codes  which  required  At  < t^.  Thus,  CP  codes  may 
be  used  to  distinguish  sharp  range  boundaries  without  any  loss  of  per- 
formaix:e. 

CP  codes  require  that  both  code  word  "A"  and  its  compliment  "B", 
illuminate  each  pixel.  Unfortunately,  the  sum  of  the  cross-correlation 
terms  of  ”A"  and  "B"  do  not  cancel.  Therefore,  care  must  be  taken  to 
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insure  that  there  is  no  cross-talk  between  the  "A"  and  "B"  transmitters 
and  receivers.  Hence  a transmitter- receiver  combination  is  required 
for  both  the  "A"  codewords  and  for  the  "B"  codewords.  The  fields-of- 
view  of  these  transmitter-receiver  pairs  must  not  overlap;  i.  e. , the  "A" 
and  "B"  channels  must  be  spatially  orthogonal.  If  not  spatially  orthogonal, 

the  "A"  and  "B"  codewords  must  be  spatially  pnd  temporally  coded  so  that 

1 

no  crosstalk  occurs  between  the  "A"  and  "B"  channels.  An  example  of 

I 

the  spatial-temporal  coding  is  given  in  the  next  chapter. 

The  CP  code  requires  two  transmitters,  and  thus  if  each  codeword 
has  as  many  chips  as  a PN  codeword,  the  received  energy  per  pixel  for  the 
CP  coding  scheme  would  always  be  twice  f^hat  of  the  PN  coding  scheme.  In 
order  to  make  comparisons  between  PN  and  CP  results,  the  CP  code  will  be 
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Fig  9.  Autocorrelation  of  a Totally  Orthogonal  Complimcntar}' 
Pair 
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required  to  have  half  as  many  chips  per  word  as  the  PN  code,  or  it  must 
have  twice  die  scan  rate  (half  the  dwell  time)  of  that  for  the  PN  code. 
Thus  the  performance  results  are  normalized  to  the  same  total  energy 
per  pixel. 


The  mean  square  bandwidth  may  be  shown  to  be 

r.a  ~ 8 


The  output  power  has  the  same  form  as  Eq  (65)  and,  therefore,  the 
detector  current  signal  energy  is  the  same  as  that  given  by  Eq  (66) 


(70) 

(71) 


where  p is  the  total  number  of  chips  in  both  A and  B.  The  mean  square 
bandwidth  and  energy  for  PN  and  CP  codes  are  approximately  the  same, 
and  thus  it  was  Eqs  (64),  (66),  and  (67)  that  were  used  to  calculate  the 
performance  for  pulse  coding  schemes  in  the  previous  section. 

Closed-Loop  Estimators.  Closed-loop  estimators  are  quite  differ- 
ent from  one-shot  matched  filters  because  they  employ  feedback.  The 
estimator  input  is  compared  to  a locally  generated  replica  of  the  signal 
and  an  error  voltage  is  generated.  The  error  signal  is  used  to  drive  the 
time  base  for  the  local  reference  in  the  appropriate  direction  to  decrease 
the  error.  Closed-loop  estimators  may  be  used  for  both  digital  and  sinu- 
soidal signals.  The  closed-loop  estimator  for  pulse  code  modulation 


would  require  many  codewords  to  slide  through  the  filter  before  the  error 
could  be  made  small.  However,  the  performance  of  the  laser  line  scan 
system  for  pulse  code  modulation  has  been  shown  to  be  best  when  the  code- 
words are  equal  to  the  dwell  time.  Therefore,  closed-loop  estimators 
would  not  be  practical  for  this  system  when  pulse  code  modulation  is  used. 

If  RF  sinusoidal  intensity  modulation  is  used,  there  would  typically 
be  many  cycles  of  the  sinusoid  in  a dwell  time.  In  this  case,  a closed- 
loop  matched  filter  is  practical.  The  best  known  closed-loop  matched 
filter  is  the  phase  locked  loop  (PLL).  If 

y(t)  = A cos(u^t  + 6 ) + n(t)  (72) 

then  it  can  be  shown  that  the  ML  estimator  for  phase  is 

y(t)  sin  (tu^t  + dt  = 0 (73) 

and  may  be  realized  by  the  PLL  below  (Fig  10)  (Ref  8:  407-409).  The 
CR  bounds  for  the  performance  of  a PLL  estimator  and  direct  detector 
are  given  by  Eqs  (55)  - (58).  From  Eq  (57)  the  range  ambiguity  for 
sinusoid  modulation  is 

^amb  “ 

Background  Noise  Considerations.  For  high  quality,  cooled 
detectors,  the  noise  rate  function  is  given  by 

" ^b  + \l 
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where 


Fig  10.  ML  Esti  rnator  for  Phase 


is  the  detected  background  noise  power  and  may  be  shown,  assuming 
a Lambertian  reflectance  target,  to  be  given  by 

Pjj  = -^  ^s.X^r"^  AX  Af  (daytime  operation)  (77) 
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where: 


Eg  = solar  spectral  irradiance  (watts /cm® -pm) 

tj.  = transmission  of  receiver  optics  and  atmospheric  path 

AX  = optical  filter  bandpass  (um) 

Aj.  = area  of  receiver  optics  (cm®) 

$ 

Af  = area  on  ground  subtended  by  receiver  FOV  (cm®) 
r = range  (cm) 

(Ref  12:  Chap  5,  p.  45).  Likewise  may  be  shown  to  be  assuming 
matched  transmitter  and  receiver  EOV's 


Pt  Aj- 

^’d  "IT  73-  ^cos  e 


(78) 


where  is  the  transmitter  power,  t^.  is  the  transmission  of  the  combined 
transmitter  and  receiver  optics  and  two  way  atmospheric  path,  and  6 is 
the  instantaneous  scan  angle  with  respect  to  the  pixel  surface  normal  vector. 
The  signal-to-noise  rate  function  ratio  is  given  by 

Xr  _ ppPd  _ Pftr  cos  9 (79) 


If  Pj  = 1 watt,  t[  /t^  « 1,  0 = 45°,  Eg  85  pm)  = 0. 007  watt/cm®  - pm 
(Ref  13:  Chap  6,  p.  2),  AX  = 0.  01  pm,  and  Af  = 2090  cm®  (18  inches  by 
18  inches),  then  Xg/Xj,  wO.  48  and  the  background  noise  limited  approxi- 
mation is  fairly  good  for  daytime  operation.  The  spectral  irradiance 
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from  the  moon  and  stars  and  the  blackbody  irradiance  of  the  earth  are 
several  orders  of  magnitude  lower  than  the  solar  spectral  Irradiance. 
Therefore,  for  nighttime  operation,  the  signal  shot  noise  limited  Poisson 
regime  is  a good  approximation. 
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III.  spatial  and  Temporal  Ccxiing  Analysis 


Introduction 

In  Chapter  II  it  was  shown  that  both  tlie  range  and  reflectance  per- 
formance variances  are  inversely  proportional  to  the  detected  signal 
energy  (Eqs  (38)  - (41)).  There  are  two  practical  ways  to  increase  tlie 
detected  signal  energy.  The  first  way  would  be  to  use  a more  powerful 
laser  source.  This  would  require  using  one  of  the  bulkier,  heavier  and 
more  fragile  solid  state  or  gas  lasers.  The  second,  and  more  attractive, 
way  would  be  to  use  more  tlian  one  GaAs  laser  diode  source.  The  rest  of 
this  chapter  will  be  devoted  to  investigating  methods  of  temporal  coding 
and  spatially  scanning  N GaAs  laser  sources  to  achieve  an  optimum  laser 
line  scan  system. 

Combining  N Source  Pignuts.  The  first  question  that  comes  to 
mind  when  discussing  the  use  of  N sources  is  "Why  not  use  a square 
n X n array  to  form  one  scanning  beam  which  will  be  detected  by  one 
scanning  receiver?",  where  n x n = N.  There  are  two  primary  reasons 
why  this  is  not  feasible.  The  first  has  to  do  witli  thermodynamics;  it  is 
very  hard  to  heat  sink  and  cool  an  array  evenly.  The  diode  lasers  at 
the  center  of  the  array  will  be  much  hotter  than  those  on  the  outside  edge 
of  the  array  and  might  be  destroyed  from  thermal  effects.  The  second 
reason  has  to  do  with  optics:  only  one  element  of  tlie  array  would  be  on 
tlie  optical  axis.  It  is  possible  to  focus  the  output  beam  of  all  of  the 
elements  onto  one  pixel  at  a given  range.  However,  as  soon  as  the 
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range  changes,  the  individual  beams  will  no  longer  be  focused  on  the 
same  pixel  and  there  will  be  interpixel  crosstalk  (noise).  A solution 
might  be  to  use  fiberoptics  to  combine  the  output  of  N laser  diodes. 

Each  diode  could  be  physically  mated  to  an  optical  fiber  which  would 
feed  into  a star  coupler.  A single  optical  fiber  out  of  the  star  coupler 
could  then  be  placed  at  the  effective  focus  of  the  system  optics.  Un- 
fortunately, the  laser  diode  to  optical  fiber  junction  has  a transmission 
loss  on  the  order  of  lOdb,  much  too  high  to  be  practical. 

The  major  part  of  this  chapter  will  be  devoted  to  investigating 
various  ways  of  using  N laser  sources  to  simultaneously  interrogate  or 
illuminate  N different  pixels  to  improve  the  overall  system  performance. 

In  order  to  gain  the  full  advantage  of  such  a multiple  source  system,  the 
N transmitter-receiver  pairs  (or  channels)  must  not  interfere  with  each 
other.  Thiais  easily  achieved  when  using  pseudo-noise  (PN)  code  modu- 
lation. Even  if  the  same  signal  energy  from  source  m falls  into  the 
field-of-view  of  detector  n,  the  amount  of  interpixel  crosstalk  will  be 
minimal  so  long  as  the  adjacent  source  modulating  signals,  Sm(t)  and 
Sjj(t),  are  very  nearly  temporally  orthogonal  (zero  cross-correlation). 

Since  the  autocorrelation  of  a PN  code  is  -1/p  everywhere  t±  o,  the 
nearly  temporally  orthogonal  condition  can  be  achieved  by  cycling  the 
PN  codeword  transmitted  by  source  n from  that  transmitted  by  source 
m by  an  amount  greater  than  the  maximum  expected  interpixel  delay 
At.  The  maximum  cross-correlation  between  Sm(t)  and  Sp(t)  will 
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then  never  be  more  than  -l/p.  Unfortunately,  complimentary  pair  (CP) 
code  modulation  is  not  orthogonal;  and  therefore,  interpixel  crosstalk 
may  not  be  minimized  by  cycling  the  codewords  of  adjacent  sources.  In 
the  CP  code  modulation  case,  interpixel  or  interchannel  isolation  must  be 
achieved  by  making  the  channel-to-channel  transmitter-receiver  pair 
fields-of-view  (FOV)  spatially  orthogonal,  i.e. , the  pixel  subtended  by 
the  FOV  of  transmitter-receiver  pair  m must  not  overlap  the  FOV  of 
transmitter-receiver  pair  n.  The  N channels  may  also  be  made 
temporally  orthogonal  by  using  a sinusoidal  modulation  of  a different 
frequency  for  each  channel.  However,  Eqs  (55),  (56)  and  (74)  reveal 
that  the  performance  and  ambiguous  range  would  vary  from  pixel  to  pixel. 
Thus,  for  sinusoidal  modulation,  a better  method  would  be  to  modulate 
all  N sources  with  the  same  frequency  sinusoid,  and  again  rely  on  spatial 
orthogonality  of  the  transmitter-receiver  FOV's  to  achieve  inter-channel 
isolation. 

Constraints  . The  performance  of  the  various  spatial -temporal 
coding  combinations  presented  in  this  chapter  will  be  expressed  using  the 
parameters  for  a one  source  line  scan  system  (Fig  H).  It  will  be  assumed 
that  the  transmitter  and  receiver  fields-of-view  are  matched,  i.  e. , subtend 
the  same  ground  area  (a  pixel)  and  are  spatially  orthogonal  from  channel 
to  channel.  The  number  of  pixels  p>er  scan  line  will  be  a fixed  integer 


number  q.  Except  for  the  compliinentary  pair  (CP)  pulse  code  modu- 
lation case,  the  dwell  time  will  be  sec  per  pixel.  In  order  to  insure 


the  same  energy  per  pixel,  the  two  sources  required  for  CP  code  modu- 
lation will  each  have  dwell  times  of  Tj/2.  The  period  of  tlie  scanline  is 
Tg  sec  where  Tg  = qT^j.  The  total  number  of  sources  will  be  N,  whether 


they  are  arranged  in  parallel,  series,  or  a m x n array  where  m x n = N. 


By  expressing  the  performance  of  each  spatial -temporal  combination  in 
these  common  terms,  performance  comparisons  can  be  easily  made  in 
Chapter  IV.  Since  Eqs  (38)  - (41)  indicate  that  the  signal  shot  noise 
regime  performance  is  improved  in  the  same  way  as  the  background  noise 
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limited  regime  performance  as  the  signal  energy  per  pixel  is  increased, 
only  the  background  noise  limited  performance  of  each  spatial-temporal 
combination  will  be  analyzed. 

N Sources  In  Parallel,  Pulse  Code  Modulation 

The  first  spatial-temporal  code  to  be  considered  places  N sources 
in  parallel  and  uses  pulse  code  temporal  modulation.  The  scanner  sweeps 
out  N parallel,  contiguous  strips  along  the  ground  (Fig  12).  This  configu- 
ration can  be  scanned  two  different  .ways.  Each  will  be  analyzed  below, 
first  using  pseudo-noise  (PN)  pulse  code  modulation,  and  then  using  comp- 
limentary pair  (CP)  pulse  code  modulation. 
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Slow  Parallel  Scan,  PN  Code  . Let  the  N sources  be  arranged 
in  a parallel  scan  configuration  and  be  temporally  modulated  by  a PN  code. 
Since  N strips  are  swept  out  instead  of  one,  the  scanner  may  be  swept  at 
1/N  the  sweep  rate  required  of  the  single  source  system  and  still  cover 
the  entire  area.  Reducing  the  scan  rate  causes  the  strips  swept  out  along 
the  ground  to  be  diagonal,  not  perpendicular,  to  the  aircraft  ground  track. 
This  should  not  pose  a serious  problem,  however,  because  any  processor 
sufficiently  sophisticated  to  provide  cueing  will  be  capable  of  de-skewing 
the  scene  raster  lines. 

Slowing  down  the  scan  by  a factor  of  1/N  is  equivalent  to  increasing 
the  dwell  time,  and  consequently  the  energy,  per  pixel  by  N times.  The 
expressions  for  the  range  and  the  reflectance  mean  square  error  for  each 
receiver  may  be  found  by  replacing  with  Nt^  in  Eqs  (50)  and  (51),  re- 
sulting in 
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c^^n^o 

C® X t 

_ ^n'-o 

(80) 

t 

N 

Spap^P^Td 

8p»p2  P^Np 

i » • 

2>^n  _ 

2^n 

(81) 

IT 

P"P§Npto 

Thus,  both  the  range  and  tlie  reflectance  mean  square  errors  are  reduced 
by  a factor  of  N.  From  Eq  (SO),  it  can  be  seen  tliat  in  order  to  realize 
the  improvement  in  range  estimation  performance,  the  pulse  codeword 
length  must  be  increased  by  a factor  of  N;  i.e. , the  nutnber  of  chips  per 
codeword  must  be  Np.  Increasing  die  chip  duration  tQ  so  that  Ntj  = pt^ 
will  not  improve  tlie  range  estimation  performaiKe. 
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The  maximum  unambiguous  range  is  proportional  to  the  codeword 
period,  from  Eq  (63),  and  becomes 


^amb 


cNTd 

“2 


(82) 


Therefore  reducing  the  sweep  of  the  scanner  also  improves  the  maximum 
unambiguous  range  by  a factor  of  N. 

Fast  Parallel  Scan,  PN  Code  . The  PN  code  temporally  modulated 
parallel  scanner  system  may  also  be  swept  at  the  same  rate  as  the  single 
source  system,  thereby  scanning  each  pixel  N times.  The  dwell  time, 
and  consequently  the  energy  per  source  per  pixel,  is  the  same  as  for  the 
one  source  system.  However,  if  the  signal  energy  received  by  each  re- 
ceiver from  a given  pixel  is  somehow  combined  before  making  the  range 
and  reflectance  estimates,  then  the  signal  energy  per  pixel  is  again  in- 
creased by  a factor  of  N.  The  N signal  returns  from  each  pixel  may  be 
combined  by  an  appropriate  network  of  matched  filters  and  delay  line 
summers  (Fig  13).  The  range  and  the  reflectance  performance  is  exactly 
the  same  as  in  the  parallel  scan  case  and  expressions  for  the  error  vari- 
ances are  given  by  Eqs  (80)  and  (81).  However,  the  dwell  time  and  the 
pulse  codeword  length  are  the  same  as  that  of  the  single  source  system, 
therefore,  the  maximum  unambiguous  range  is  given  by 
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Fig  13.  Receiver  Configuration  for  "Fast  Parallel 
Scan  - PN  Code"  System 


Besides  the  shorter  unambiguous  range,  this  configuration  requires 
q of  the  relatively  long  delay  line  summers  shown  in  Figure  11,  one  for 
each  pixel  in  the  scan  line.  However,  it  does  have  the  advantage  that 
any  differences  in  individual  laser  power  outputs  and  individual  detector 
responsivities  are  averaged  out  by  the  delay  line  summers.  Also,  the 
strips  swept  out  on  the  ground  will  be  much  closer  to  perpendicular  to 
the  aircraft  ground  track  than  those  in  the  parallel  scan  case. 

Slow  Parallel  Scan,  CP  Code  . Now  consider  the  N parallel  scan- 
ning sources  temporally  modulated  with  a CP  code.  Let  the  first  N/2 
sources  transmit  codeword  "A",  the  second  N/2  sources  transmit  the 
complimentary  codeword  "B",  and  the  N sources  be  scanned  at  a rate 
which  is  2/N  that  of  the  one  source  system  scan  rate.  Each  pixel  will 
be  scanned  twice,  first  by  a source  transmitting  codeword  "A",  and  then 
NTg/2  sec  later  by  a source  transmitting  the  complimentary  codeword 
"B",  where  T„  is  the  scan  period  of  the  single  source  system.  Now,  the 

o 

effective  dwell  time  is  NTjj/2  and  each  pixel  is  scanned  twice  so  that  the 
total  received  signal  energy  per  pixel  is  again  increased  by  a factor  of  N. 
The  two  signal  returns  from  each  pixel  may  be  processed  by  using  a delay 
line  (Fig  14).  The  performance  is  exactly  the  same  as  in  the  previous 
two  cases  and  the  expressions  for  the  error  variances  are  given  by  Eqs 
(80)  and  (81).  From  Eq  (63)  the  maximum  unambiguous  range  is 
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i.  e. , N/2  times  that  for  the  single  source  system. 


This  spatial-temporal  combination  requires  qN/2  delays,  has  a 
skewed  scene  raster,  and  only  averages  responses  of  two  transmitter- 
detector  pairs.  It  also  requires  strict  spatial  orthogonality. 


Fig  14.  Receiver  Configuration  for  CP  Coding 


Fast  Parallel  Scan,  CP  Code  . The  N sources  may  also  be  scanned 
at  the  same  rate  as  the  single  source  system  while  being  temporally 
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modulated  by  a CP  code.  The  only  difference  between  this  and  the  PN 
case  discussed  above  is  that  one-half  of  the  sources  must  transmit  code- 
word "A"  and  the  other  half  must  transmit  codeword  "B”.  There  are 
several  ways  this  may  be  done,  but  one  way  might  be  to  let  all  of  the  odd 
numbered  sources  transmit  codeword  "A"  and  let  all  of  the  even  numbered 
sources  transmit  codeword  "B".  The  receiver  outputs  may  be  combined 
by  using  the  configuration  shown  in  Figure  11  where  all  of  the  odd  matched 
filters  are  matched  to  codeword  "A"  and  all  of  the  even  matched  filters 
are  matched  to  codeword  "B”.  The  performance  expressions  and  unam- 
biguous range  are  exactly  the  same  as  for  the  "Fast  Parallel  Scan-PN  Code" 
case  and  are  given  by  Eqs  (80),  (81),  and  (83).  The  advantages  and  dis- 
advantages are  nearly  identical  also.  This  case  has  the  additional  dis- 
advantage of  requiring  N more  delays,  and  has  the  additional  advantage 
that  the  CP  code  is  linear  and  will  not  be  affected  by  inter-pixel  delay 
differences. 


N Sources  In  Series,  Pulse  Code  Modulation 

Another  possible  spatial-temporal  code  places  N sources  in 
series  and  uses  pulse  code  temporal  modulation  (Fig  15).  Since  only 
one  strip  of  ground  is  swept  at  a time,  the  scan  rate  must  be  the  same 
as  that  of  the  single  source  system  in  order  for  the  scanned  strips  to  be 
contiguous.  Thus,  as  in  the  "Fast  Parallel  Scan"  case,  the  energy  per 
source  per  pixel  is  the  same  as  the  single  source  case;  therefore,  once 
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again,  all  of  the  receiver  outputs  must  be  combined  before  making  range 
or  reflectance  estimates.  If  CP  code  modulation  is  used,  with  alternating 
codewords  "A"  and  "B",  the  N receiver  outputs  may  be  combined  via  the 


Fig  15.  N Sources  in  Series 


receiver  configuration  shown  below  (Fig  16).  If  PN  code  modulation  is 
used,  the  matched  filters  are  merely  matched  to  the  PN  codeword.  It 
can  easily  be  shown  that,  just  as  in  the  parallel  scan  case,  the  performance 
is  the  same  regardless  of  whether  PN  or  CP  code  modulation  is  used,  and 
is  equal  to  that  of  the  parallel  scanning  systems.  Also,  since  the  dwell 
time  is  for  the  series  scan  regardless  of  the  type  modulation  used,  the 
maximum  unambiguous  range  is  also  the  same  as  that  of  the  single  source 

I 
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system.  Thus,  the  performance  expressions  and  unambiguous  range  are 
given  by  Eqs  (80),  (81),  and  (83).  This  system  has  one  large  advantage 


over  the  parallel  scan  systems.  In  addition  to  having  nearly  perpendicular 
raster  lines  and  averaging  out  the  responses  of  the  different  transmitter- 
detector  pairs,  only  one  relatively  short  delay  line  summer  is  required  to 
combine  the  N detector  outputs. 


N Sources  In  An  Array,  Pulse  Code  Modulation 

It  is  conceivable  that  a N source  scanner  might  be  configured 
using  a combination  of  the  parallel  and  series  scanning  systems  and 
temporally  modulated  with  pulse  code  modulation.  It  can  easily  be 
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shown  that,  so  long  as  the  total  number  of  sources  is  N and  all  of  the 
detector  outputs  are  summed  before  processing,  then  the  performance 
of  this  type  system  will  be  exactly  that  of  either  the  parallel  or  series 
3can  systems,  and  will  be  independent  of  which  type  of  code,  PN  or  CP, 

Is  used. 

The  major  things  that  will  vary  with  configuration  will  be  the  maxi- 
mum unambiguous  range,  the  complexity  of  the  receiver  configuration,  the 
angle  of  the  scene  raster  lines  with  respect  to  the  ground  track,  and  the 
number  of  transmitter-detector  pair  responses  that  are  averaged.  For 
an  m X n array,  let  there  be  m sources  in  parallel,  and  n sources  in 
series.  From  the  discussions  above,  it  is  clear  that  the  simplest  delay 
line  summer  network  will  be  realized  if  the  system  is  scanned  at  the 
"slow"  rate  of  1/m  tliat  of  the  single  source  system,  and  no  detector  out- 
puts are  combined.  In  tins  configuration,  a block  of  m contiguous  strips 
are  swept  out,  and  then  the  system  moves  forward  and  sweeps  out  another 
block  of  m strips  adjacent  to  the  first  block.  Once  more  the  range  and  ■ 

reflectance  performance  expressions  are  given  by  Eqs  (80)  and  (81). 

The  dwell  time  is  mi^,  and  thus  the  maximum  unambiguous  range  is 

p _ cmT(j  (85) 

^amb  — ^ — 

For  this  configuration,  m delay  line  summers  are  required  to  combine 
the  n detector  outputs  per  scan  line.  The  greatest  disadvantage  of  this 


52 


configuration  is  that  an  m x n array  of  laser  diodes,  either  on  one  sub- 
strate or  in  an  array  of  individual  diodes,  is  required  for  the  transmitter. 
The  former  is  unlikely  because  of  the  difficulty  in  cooling  the  array;  and 
the  latter  is  undesirable  due  to  the  bulk  and  complicated  optics  involved. 

N Sources  In  Parallel,  Sinusoidal  Modulation 

The  N parallel  source  scan  configuration  shown  in  Fig  10  may  also 
be  sinusoidally  temporally  modulated.  If  the  scanner  is  swept  at  1/N  the 
rate  of  the  single  source  system,  then  die  dwell  time  is  once  again  Nt^. 
The  expressions  for  the  range  and  the  reflectance  mean  square  errors  may 
be  found  by  dividing  Eqs  (55)  - (56)  by  N,  resulting  in 

^ JL  • - (86) 
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where  f^^  is  the  frequency  of  die  sinusoidal  intensity  modulation  of  the 
N sources. 

The  maximum  unambig^uous  range  is  given  by  Eq  (74),  and  is 


^amb  ~ 2frn 


Note  diat  the  performance  is  inversely  proportional  to  fj^  squared  where 
as  the  range  ambiguity,  which  is  independent  of  dwell  time  for  sinusoidal 
modulation,  is  inversely  proportional  to  f_.  Thus,  there  is  a definite 
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tradeoff  between  system  performance  and  maximum  unambiguous  range. 

For  example,  if  the  laser  were  intensity  modulated  near  its  maximum 
frequency  response,  say  lOOMlIz,  to  minimize  the  range  and  reflectance 
errors,  the  unambiguous  range  would  be  1.  5 meters  (4.  9 feet). 

This  coding  combination  will  also  exhibit  the  skewed  raster  lines 
as  in  the  "Slow  Parallel  Scan,  Pulse  Code"  cases  discussed  earlier. 

Another  disadvantage  is  that  any  differences  in  transmitter-detector  pair 
responsivity  will  not  be  averaged  out  as  in  the  fast  scanning  configuration. 

For  an  implementation  of  a fast  parallel  scan  configuration  or  a 
series  scan  configuration  using  sinusoidal  modulation,  it  can  be  shown  that 
the  performance  and  maximum  unambiguous  range  is  just  that  of  the  slow 
parallel  scan  case.  The  advantages  of  the  fast  parallel  or  serial  scans 
would  be  the  reduced  skewing  of  the  raster  lines  and  the  averaging  of  the 
N channel  responsivities.  Unfortunately,  these  implementations  are  not 
practical.  As  in  the  pulse  code  modulation  cases  above,  it  would  be  neces- 
sary to  coherently  corT±)ine  the  N receiver  outputs  from  each  pixel  before 
making  range  or  reflectance  estimates.  The  coherent  addition  of  RF 
frequencies  using  many  fairly  long  delay  lines  would  be  very  difficult. 

N Sources  In  Parallel,  Multiple  Sinusoidal  Modulation 

Using  existing  technology,  sinusoidal  in  tensity  modulation  of  high 
powered  GaAs  laser  diodes  can  be  done  (Ref  2).  For  this  reason,  it  is 
highly  desirable  to  find  some  kind  of  sinusoidal  temporal  modulation 
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which  has  a large  maximum  unambiguous  range.  One  such  method  is  to 
use  two  superimposed  sinusoids  of  different,  but  close,  frequencies 
(Fig  17). 
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Fig  17.  Multiple  Sinusoid  Transmitter  Configuration 
If  the  two  sinusoids  are  given  by 

Si(t)  = sin(2nf^t+$p  (88) 


S^(t)  = sin  (2TTfgt+ $3) 


then  the  range  and  maximum  unambiguous  range  are 
c A i 
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where  A $ is  the  phase  difference  between  the  two  received  sinusoids  and 

f = f + M (Ref  14:  106-107).  Thus  such  a modulation  scheme  has  a 
a 1 

potential  for  coarse  range  measurement  with  long  unambiguous  range 
provided  that  the  standard  deviation  of  Af  is  less  than  Fine  range 
resolution  nrea  sure  men  ts  can  be  nnade  by  routing  the  output  of  a high-pass 
filter  centered  at  f^  into  a phase  locked  loop. 

The  performance  of  the  fine  range  estimates  may  be  found  from 
Eq  (55),  which  may  be  expressed  as 

oi  _£lis (91) 

16  TT^f*  p^Eg 

where  Eg  is  given  by  Eq  (48).  Since  the  laser  diodes  are  assumed  to  be 
peak  power  limited,  the  vector  sum  of  the  two  modulating  sinusoids  may 
not  exceed  the  laser's  peak  output  power.  Therefore,  s(t)  in  Eq  (54) 
becomes 

PPd 

s(t)  = —2 — (1  + 1 sin  2TTfit  + h sin  2TTfat)  (92) 

Substituting  the  high-pass  filtered  output,  (pp  P^/4)  sin  2TTf2t,  and  Eq 
(48)  into  Eq  (90)  yields 

i (93) 

The  performance  of  the  reflectance  estimates  may  be  found  by 
substituting  Eqs  (92)  and  (48)  directly  into  Eq  (39)  yielding 
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Dividing  Eqs  (93)  - (94)  by  N and  then  comparing  to  Eqs  (86)-(87) 
reveals  that  the  effect  of  splitting  the  available  laser  power  between  two 
sinusoids  to  increase  range  ambiguity  by  a factor  of  fg/Af  is  to  increase 
the  mean  square  range  error  by  a factor  of  (4)  and  the  mean  square  reflect- 
ance error  by  a factor  of  (1.  2). 

Temporal  Modulation  Via  Electro-optic  Modulators 

High  power  laser  diodes  have  low  input  impedances  and  require  high 
drive  currents.  Thus,  narrow  pulses  are  difficult  to  obtain  because  the 
extremely  large  voltage  drops  across  associated  circuit  inductances  as 
small  as  a nano-henry.  Since  Eq  (71)  implies  that  the  range  performance 
of  pulse  code  modulation  schemes  is  proportional  to  the  chip  width  t^, 
indirect  methods  of  obtaining  narrow  intensity  pulses  are  attractive. 
Temporal  intensity  modulation  may  be  achieved  using  electro-optic  modu- 
lators such  as  Kerr  cells  or  Pockel  cells  (Ref  15:  245-268).  Gigabit  per 
second  modulation  rates  with  24  volt  modulation  voltage  and  81.  5%  static 
transmission  for  530nm  light  have  been  obtained  using  Pockel  cells 
(Ref  16). 

If  the  high  power  GaAs  laser  diodes  are  assumed  to  be  unpolarized, 
then  tliere  will  be  an  automatic  3 db  power  reduction  due  to  the  Pockel  cell 
polarization  analyzer.  Thus,  the  total  through-put  would  be  approximately 
0.4  Pj.  Assuming  peak  power  limited  laser  diodes,  the  previous  pulse 
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modulation  analysis  may  be  used  by  replacing  everywhere  by  0.4  P^. 

I Thus  the  range  and  reflectance  rms  errors  will  be  2. 5 times  larger  than 

I;  the  case  of  identical  parameters  obtained  by  direct  modulation  of  the 

I'  laser  diode  drive  current.  However,  if  the  electro-optic  modulator  is 

capable  of  chip  widths  less  than  0. 16  that  available  by  direct  modulation, 
then  the  loss  of  power  will  have  been  overcome,  and  the  electro-optic 
modulator  will  provide  the  best  performance.  This  result  is  valid 
assuming  that  the  overall  system  response  is  limited  by  the  speed  of  the 
j transmitter,  and  not  that  of  the  detector. 

' Temporal  Modulation  Via  Spatial  Coding 

Another  indirect  way  of  generating  a pulse  at  the  output  of  the 
1 receiver  is  to  stare  at  a given  pixel  with  a detector  and  sweep  the  trans- 

(i  mitted  beam  from  a CW  source  through  the  detector  field-of-view  (FOV). 

j Thus,  it  is  possible  to  temporally  modulate  the  energy  return  from  a pixel 

I by  spatially  codi  ng  the  scanner.  The  detector  output  is  just  the  spatial- 

I temporal  convolution  of  the  transmitter  and  receiver  FOVs.  If  the  scan- 

; ning  transmitter  beam  is  spatially  coded  with  a given  sequence,  then  the 

1 detector’s  output  will  look  like  a pulse  code  word  with  the  same  sequence. 

This  may  be  accomplished  two  ways:  (1)  N-series,  spatially  coded 
sources  may  be  scanned  p/2  times  per  scan  line  by  a rotating  polygonal 
prism  with  p/2  faces  (Fig  18);  and  (2)  N-CW  sources  in  parallel  may  be 
scanned  by  a coded  rotating  polygonal  prism  with  p faces: 
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Fig  18.  N Series,  Spatially  Coded  Sources 

a digital  "one"  would  correspond  to  a highly  polished  prism  face,  and  a 
digital  "zero"  would  correspond  to  a rough  prism  face  (Fig  19).  The 
length  of  the  codeword  is  a 2N  chips  for  the  first  case  and  p chips  for 
the  second  case.  So  long  as  the  number  of  sources,  the  number  of 
sides  of  the  polygon,  and  the  chip  duration  (dwell  time)  t^^  are  such  that 
NTj  = NptQ  then  the  rms  errors  for  the  range  and  reflectance  estimates 
will  be  given,  once  again,  by  Eqs  (80)  and  (81),  and  the  unambiguous 
range  will  be  given  by 

^amb  ” 
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Fig  19.  N Parallel,  Spatially  Coded  Transmitter  Beams 


for  N series,  spatially  coded  sources,  and  by 


R 


cpto 


amb 


(96) 


for  N sources  in  parallel.  Hence,  if  t^^  can  be  made  shorter  via  spatial 
coding,  then  the  performance  will  be  improved.  The  price  that  is  re- 
quired for  this  improvement  is  the  need  for  ultra-fast  scanner  rotation 
rates,  large  numbers  of  detector  elements,  and  complicated  optics.  It 
should  be  noted  that  a beam  splitter  would  be  required  to  sample  the  trans- 
mitted beam  in  order  to  obtain  a reference  from  which  to  measure  delay. 

A very  small  amount  of  optical  signal  (from  the  beam  splitter)  would  be 
detected  and  processed  by  the  matched  filter.  The  time  between  the 


60 


occurrence  of  a peak  at  the  matched  filter  output  due  to  the  reference  and 
the  occurrence  of  a peak  due  to  the  reflected  signal  is  the  delay  estimate 


Hybrid  Spatial  and  Temporal  Coding 


In  an  effort  to  increase  the  bandwidth  of  the  intensity  modulation,  a 
transmitted  laser  beam  may  be  both  temporally  modulated  and  spatially 
coded,  thereby  producing  a hybrid  modulation  at  the  detector  output.  If 
pulse  code  temporal  modulation  is  used  and  the  codeword  period  is  longer 
than  a dwell  time,  the  effect  will  be  to  modulate  or  "gate"  the  temporal 
modulation  with  the  spatial  code  (Fig  20).  Thus,  if  the  autocorrelation 
properties  of  the  PN  or  CP  codes  are  to  be  preserved,  the  codeword  dura- 
tion must  be  less  than  or  equal  to  the  dwell  time.  Hence  there  is  no  ad- 
vantage to  gating  a pulse  code  temporal  modulation  via  spatial  coding. 


Fig  20.  Hybrid  Spatial-Temporal  Code 
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However,  if  the  source  is  temporally  modulated  with  a sinusoid, 
the  result  of  spatial  coding  will  be  to  On-Off-Key  (OOK)  the  sinusoid.  If 
the  spatial  code  is  a pulse  code,  then  the  output  of  the  photodetector  is  a 
pulse  coherent  pulse  code  (Fig  21).  Thus,  if  a phase  locked  loop  (PLL) 
can  achieve  phase  lock  by  the  end  of  the  pulse  codeword,  it  may  then  be 
used  to  make  range  and  reflectance  estimates:  and  the  pulse  codeword  may 
be  coherently  detected  and  used  to  resolve  the  range  ambiguity  of  the  sinu- 
soid. 

Recall  that  the  reason  for  considering  the  spatial  coding  in  the  first 
place  was  to  obtain  a way  of  creating  pulses  with  pulsewidths  smaller  than 
could  be  achieved  via  direct  modulation  the  laser  diode  drive  current. 


mNwmiwmNwmim 

Sinusoidal  Temporal  Modulation 

PN  Code  Spatial  Modulation 

,A/VVVWVVWAVv'MW 

Photodetector  Output 

WVvAM 

Fig  21.  On-Off-Keying  Via  Hybrid  Spatial-Temporal  Coding 

The  OOK  sinusoidal  modulation  does  not  require  narrow  chips.  In  fact, 

the  chips  must  be  long  enough  for  many  cycles  of  the  sinusoid  to  fit  inside  i 

j 

I 
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a chip  so  that  the  PLL  can  achieve  phase  lock.  Hence,  there  is  no  advan- 
tage to  gating  sinusoidal  temporal  modulation  via  spatial  coding.  The 
disadvantages  include  high  scanner  rotation  rates,  large  numbers  of  detect- 
or elements,  and  complicated  optics. 

N Sources  in  Parallel,  OQK  Sinusoidal  Modulation 

Each  source  of  the  N parallel  source  scan  configuration  shown  in 
Fig  12  may  be  CXDK  sinusoidal  temporally  modulated  by  direct  modulation 
of  the  laser  source  drive  current  (Fig  22).  The  range  and  reflectance  esti- 
mates for  each  of  the  N detectors  may  be  made  via  the  receiver  configuration 


below  (Fig  23).  The  voltage  controlled  oscillator  (VCO)  constant 
(see  Fig  16)  must  be  small  in  order  for  the  PLL  to  acquire  and  maintain 
phase  lock  with  OOK  sinusoidal  input.  Thus,  two  PLLs  are  alternated 
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Reset 


Fig  23.  Receiver  Configuration  for  OOK  Modulation 


to  make  the  phase  estimates.  A local  oscillator  samples  and  "remembers" 
the  phase  and  frequency  of  the  PLL  at  the  end  of  a dwell  time  (codeword 
period).  The  local  oscillator  output  is  then  used  to  coherently  detect  the 
delayed  OOK  pulse  codeword.  Each  PLL  is  reset  before  being  used  to 
make  phase  estimates  for  a particular  pixel. 

The  Cramer  Rao  performance  bounds  for  the  mean  square  range 
and  reflectance  errors  may  be  found  from  Eqs  (86)  and  (87).  The  OOK 
with  a pulse  code  has  the  effect  of  cutting  the  illumination  time  of  each 
pixel  in  half,  thereby  reducing  the  signal  energy  out  of  the  detector  by  one- 
half;  thus,  replacing  T^by  r^/2  in  Eqs  (86)  and  (87), 


1 . c * ^n 

N n^f  2p 


(97) 
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(98) 


The  resulting  mean  square  errors  for  range  and  reflectance  are  2 times 
larger  than  for  sinusoidal  modulated  sources. 

The  input  to  the  pulse  code  matched  filter  is 


Low  Pass  { p p P(t)  [1  + cos  ( “'m*^  + 0 ) } 

= P (99) 

where  P(t)  is  given  by  Eq  (65)  with  Pj.  replaced  by  P^j . 
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Therefore,  the  signal  energy  available  to  tlie  matched  filter  is  just  one 
fourth  wliat  it  was  for  the  pulse  code  modulation  case;  thus 


a 


a 

l^amb 


c=Xnto 
2 p3  p 


(100) 


A conventional,  though  a rbitinry,  measure  of  range  resolution  is 
that  two  target  echoes  must  be  separated  by  at  least  tlie  full  width  at  half 
maximum  of  the  matched  filter  output.  Similarly,  in  order  for  die  pulse 
code  to  help  resolve  die  sinusoidal  range  ambiguity,  a logical,  though  arbi- 
trary, requirement  is  that 


2 Od  ^ ^ -rjT — 

‘'a  mb  -^‘m 


(101) 


Substituting  Eq  (100)  into  Eq  (101)  and  solving  for  t^  yields 


® ^n^m 


(102) 


For  example,  if  the  sinusoidal  modulation  frequency  is  100  MHz,  = 

0.  1 X„  (background  noise  limited  case),  Pjj  = 10"®  watts,  p = 2.  5 x 10^® 
jouk‘8"'  and  p = 1.3,  then  Iq  must  be  less  than  or  equal  to  41  nscc,  approxi- 
rnarrlv  f<Hir  tiiiK’s  tiic  perioil  of  die  sinusoid.  As  die  signal- to- noise  ratio 
•«vi  •»f  '♦!#•  miMil'i'i  »•(  chips  p ini.re.i.ses.  die  longer  t^  maybe  while  sdll 
»-  *•>  ■•»  -I*  Ml  •■vprt-siN’d  h\  I'q  (102).  ('Iiips  widths  of  this 
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laser  dicxte  and  associated  drive  circuitry  has  a natural  resonant  frequency 
near  , which  is  usually  the  case. 

The  maximum  unambiguous  range  for  the  OOK  sinusoidal  modula- 
tion is  the  same  as  for  the  N Sources  in  Parallel,  Slow  Scan,  PCM  case, 
and  is  given  by  Eq  (82) 

R K = N £!<i  t=  Ncptp 

amb  2 (82) 

It  is  not  practical  to  try  to  implement  the  OOK  sinusoidal  modula- 
tion with  a fast,  parallel  source  scan  configuration,  or  with  a serial  source 
scan  configuration.  This  is  because  the  signal  energy  must  be  above  a 
certain  threshold  fora  PLL  to  acquire  and  maintain  phase  lock.  Thus, 
the  output  OOK  sinusoidal  modulation  from  tlie  N individual  receivers  for 
each  pixel  must  be  coherently  combined  before  being  applied  to  the  PLL. 

As  stated  earlier,  the  coherent  addition  of  many  RF  frequency  signals  would 
be  extremely  difficult. 

q- Staring  Sources  and  [\?tcctors,  Pulse  Code  Modulation 

A non-scanning  sensor  configuration  may  be  realized  with  q source- 
detector  pairs  (one  for  each  pixel  of  the  scan  liixj).  Although  this  scheme 
requires  more  sources  than  the  previous  cases,  it  has  the  advantage  of  no 
moving  parts:  tlie  staring  source-detector  pairs  swept  out  tlie  scene  as  the 
aircraft  moves  forward  (Fig  24).  The  effective  dwell  time  is  just  the  time 
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it  takes  for  the  aircraft  to  move  forward  a pixel  width,  and  is  q times  the 
dwell  time  for  the  single  source,  The  length  of  the  pulse  codeword 

may  now  be  quite  long  to  provide  very  good  estimation  performance.  For 
pulse  code  modulation,  either  PN  codes  with  length  qr^  or  CP  codes  with 
two  words  of  length  q'rjj/2  , the  expression  for  the  range  and  the  reflectance 
mean  square  error  become 


Fig  24.  q - Staring  Source- Detector  Pairs 


The  maximum  unambiguous  range  becomes 


for  PN  code  modulation,  and 


^amb-CP  " ^ “4” 

for  complimentary  code  modulation.  Comparing  these  results  with  those 
given  by  Eqs  (71)  - (73)  and  (75)  for  N pulse  code  modulated  sources,  the 
range  and  the  reflectance  mean  square  errors  and  the  maximum  unambigu- 
ous range  are  all  improved  by  a factor  of  q/n.  The  price  required  for  this 
improvement  is  q/n  times  more  source- detector  pairs  and  the  inability  to 
average  out  the  differences  in  source -detector  responsivities. 

q-  Staring  Sources  and  Detectors,  Spatial  and  Temporal  Modulation 

The  q-staring  sources  described  in  the  previous  section  may  be 
spatially  coded,  or  may  be  scanned  conventionally,  by  electronic  means. 

In  order  to  separate  the  return  from  each  pixel,  it  is  required  that  the  q 
source- detector  pairs  be  spatially  or  temporally  orthogonal,  or  both. 

Spatial  orthogonality  maybe  hard  to  achieve  using  GaAs  lasers  for  sources 
because  they  have  rather  large  beam  divergences.  Thus,  in  order  to  re- 
duce channel-to-channel  cross-talk  as  much  as  possible,  it  may  become 
necessary  to  restrict  the  source  modulation  so  that  no  adjacent  sources 
transmit  simultaneously.  Such  a scheme  is  described  below  for  CP  code 
modulation. 

Suppose  the  projection  of  the  beam  is  made  smaller  tlian  the  pixel, 
and  each  pixel  is  divided  into  five  intervals  and  die  q-elemcnt  linear  array 
is  divided  into  5 different  element  groups  (Fig  25).  Each  element  of  a 
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given  group,  group  1 (all  of  the  sources  labeled  "1"  in  Fig  25)  for  instance, 
are  separated  by  four  elements  from  the  other  four  groups.  All  of  the 
elements  of  group  1 simultaneously  transmit  codeword  "A"  of  the  CP  code, 
and  then  transmit  codeword  "B"  of  the  CP  code.  Then  all  of  the  elements 


of  group  2 simultaneously  transmit  codeword  "A"  and  then  codeword  "B", 
and  so  on.  The  effect  is  to  produce  the  spatial-temporal  ground  coverage 


shown  below  (Fig  25).  The  effective  dwell  time  for  each  pixel  is  now 


Fig  25.  q-Staring  Source- Detector  Pairs,  Spatially 
and  Temporally  Modulated 
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2qT^/5 , and  the  range  and  the  reflectance  mean  square  errors  are  given 
by 
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(107) 
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and  the  maximum  unambiguous  range 


^amb 


qci^ 
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is  given  by 


(109) 


This  particular  method  of  spatial  coding  into  the  q source-detector  pairs 
results  in  2.5  times  degradation  of  both  performance  bounds  and  the 
ambiguous  range  when  compared  to  the  previous  coding  scheme. 

Odier  configurations  with  more  than  five  different  element  groups 
may  be  arrived  at.  Also  if  peak  power  limited  output  restraint  is  dropped 
and  rms  output  is  the  same  regardless  of  the  duty  cycle  - then  performance 
is  the  same  as  the  previous  case. 


N Overlapping  Sources , Pulse  Code  Modulation 

The  final  spatial-temporal  coding  scheme  to  be  discussed  in  the 
chapter  considers  the  spatial  combination  of  N sources  by  overlapping  the 
N transmitted  laser  beams  (Fig  26).  No  attempt  is  made  to  tightly  colli- 
mate the  beams.  Instead,  they  are  allowed  to  diverge  rather  broadly  to 
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form  a large  "spotlight"  type  of  "footprint"  on  the  ground.  Except  at  the 
edges,  the  average  radiant  emittance  (watt/m^)  per  pixel  for  such  a spatial 
distribution  is  N times  that  for  a single  source  line  scan  system.  An  array 
of  N spatially  orthogonal  detectors  are  used  to  partition  the  large  "foot- 
print" into  an  array  of  N pixels  (Fig  27).  The  "footprint"  and  N detectors 
FOVs  may  be  scanned  together,  and  all  N of  the  sources  may  be  temporally 
modulated  simultaneously  with  the  same  pulse  code.  So  long  as  the  output 
of  the  detectors  are  combined  so  that  the  differences  in  radiant  emittance 
per  pixel  and  detector  responsivities  are  averaged  out,  and  not  too  much 
energy  overspills  the  FOV  of  the  N receivers,  the  performance  can  be 
shown  to  be  exactly  that  of  the  other  spatial-temporal  coding  combinations 
which  used  pulse  code  modulation.  The  slow  scan  mxn  array  receiver 
configuration  is  die  simplest  to  implement.  However,  a fast  scanning 
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Fig  27.  "Footprint"  of  N-Overlapping  Beams  Partitioned 
By  N“ Detector  FOVs. 

mxn  array  receiver  which  scans  each  line  m times  and  uses  delay  lines  to 
sum  the  mxn  returns  from  each  pixel  will  give  the  most  uniform  reflect- 
ance performance.  The  expressions  for  the  range  and  the  reflectance  vari- 
ances are  given  by  Eqs  (80)  and  (72)  and  the  unambiguous  range  is  given  by 


(110) 


where  is  the  dwell  time  (codeword  length)  of  each  detector.  The  advan- 
tage of  such  an  illumination  scheme  is  that  the  requirements  on  the  trans- 
mitter optics  (usually  diffraction  limited)  are  greatly  relaxed,  and  the 
differences  in  all  the  transmitter-powers  and  m of  the  photodetector  res- 
ponsivities  are  averaged  out.  Also,  the  number  of  transmitting  laser 
diodes  is  arbitrary  and  need  not  be  matched  to  the  number  of  photodetectors 
used  to  partition  the  receiver  FOV.  The  number  of  sources  could,  theor- 
etically, be  several  times  that  of  the  detectors.  However,  in  order  to 
make  systems  comparisons  easier,  the  number  of  sources  will  be  limited 
to  N for  this  analysis. 
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IV.  Comparisons  of  the  Various  Spatial  -Temporal 
Coding  Combinations 

The  performance  of  the  various  spatial -temporal  coding  combi- 
nations developed  in  Chapter  III  will  be  compared  in  this  chapter.  The 
basis  of  comparison  to  be  used  will  be  that  of  the  range  and  the  reflectance 
estimation  performance  of  the  coding  schemes.  Comparisons  will  be 
made  by  relating  the  rms  (standard  deviation)  of  the  range  and  the  reflect- 
ance errors  for  each  scheme.  In  addition,  the  maximum  unambig  <us 
range  and  the  ease  of  implementation  of  the  systems  will  be  compared. 
Throughout  all  of  the  comparisons,  it  will  be  assumed  that  the  values  of 
the  laser  output  power  Pj.  and  the  photodetector  responsivity  P = “n/hfQ 
for  each  of  the  N laser  transmitter-photodetector  pairs  are  identical  from 
pair-to-pair  and  coding-scheme  to  coding-scheme.  It  will  also  be  assumed 
that  the  GaAs  laser  diodes  are  peak  power  limited  and  that  the  system  modu- 
lation bandwidth  is  limited  by  the  high  frequency  response  of  the  laser  diode 
or  the  phot  odetector,  whichever  is  smaller. 

In  order  to  simplify  the  comparisons,  the  various  spatial-temporal 
combinations  presented  in  Chapter  III  will  be  divided  into  the  three  following 
classes:  (1)  Scanning  configurations  using  pulse  code  modulation  (PCM), 

(2)  Scanning  configurations  using  sinusoidal  modulations,  and  (3)  Non- 
scanning or  staring  configurations.  Comparisons  will  be  made  within  each 
class  to  determine  the  coding  scheme  with  the  best  performance.  The 
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resulting  systems,  one  or  two  from  each  class,  will  then  be  compared. 

Also,  since  it  was  shown  in  Chapter  III  that  there  is  no  advantage  to  em- 
ploying "Hybrid  Spatial  and  Temporal  Coding",  it  will  not  be  included  in 
the  following  comparisons. 

It  should  be  noted  that  the  range  performance  will  be  the  primary 
criterion  for  determining  one  system  over  another.  As  was  shown  by  Eq 
(39),  the  reflectance  performance  depends  only  upon  the  postdetection  sig- 
nal energy  and  will  be  identical  for  all  systems  given  the  same  total  energy 
transmitted  per  pixel.  It  can  also  be  shown  that  the  range  performance  of 
a system  will  degrade  faster  than  the  reflectance  performance  as  the  signal- 
to-noise  ratio  (SNR)  is  decreased.  The  range  estimation  is  made  using 
either  a matched  filter  of  a phase  locked  loop;  and  both  require  the  signal- 
to-noise  ratio  to  be  above  some  threshold.  However,  the  reflectance 
estimation  is  made  by  simply  integrating  the  receiver  output,  and  is  possible 
at  SNRs  much  lower  than  the  threshold  required  for  operation  of  the  matched 
filter  or  PLL.  Also,  the  reflectance  data  is  of  secondary  importance  to  the 
pattern  recognition  algorithms  which  primarily  require  range  contours 
(derived  from  range  gating). 

Comparisons  of  the  Pulse  Code  Modulation  Class 

All  of  the  scanning  configurations  that  use  pulse  code  modulation, 

"N  Sources-Parallel/Series/Array/Overlapping  Beams",  have  identical 
expressions  for  the  range  and  the  reflectance  variances.  Thus  the  rms 


75 


error  expressions  will  also  be  identical.  This  is  true,  regardless  of 
whether  the  PCM  modulation  is  the  result  of  direct  modulation  (via  direct 
current  modulation)  or  indirect  modulation  (via  electro-optic  or  spatial 
modulation).  The  major  differences  between  these  scan  configurations 
are  the  resulting  maximum  unambiguous  range,  the  implementation  com- 
plexities, and  the  number  of  averaged  transmitter- receiver  pair  responsi- 
vities. 

Maximum  Unambiguous  Range.  The  maximum  unambiguous  range 
depends  upon  the  codeword  length,  or  pixel  dwell  time,  and  thus  is  a func- 
tion of  the  scan  configuration.  Comparing  Eqs  (82)-(83)  yields 

^amb-Slow  Parallel  Scan,  PN  Code 

~ ^amb-Single  Source  Scanner,  PCM 
^amb-Fast  Parallel  Scan,  PCM 
Parallel  Scan,  CP 

^amb-Scanning  Array,  PCM 
m 

— ^a mb -Overlapping  Beams,  PCM 
m 

Thus  the  slow  parallel  scan  is  capable  of  the  largest  unambiguous  range 
estimate.  However,  it  should  be  noted  that  for  a single  source  scanner 
with  Tjj  = liisec,  the  unambiguous  range  is  150  meters  (490  feet),  which 
should  be  adequate  for  providing  target  cues  with  a downward  looking  sys- 
tem. Therefore,  all  of  the  systems  using  pulse  code  modulation  should 
provide  adequate  maximum  unambiguous  ranges  for  downward  looking 
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systems.  If  the  laser  line  scanner  is  to  be  used  in  a forward  looking  con- 
figuration where  long  slant  ranges  are  involved,  then  longer  unambiguous 
ranges  may  be  desirable  and  the  slow  parallel  or  array  scanning  configu- 
ration might  be  more  desirable.  Otherwise,  maximum  unambiguous  range 
is  not  a key  parameter  to  be  used  in  recommending  one  of  the  above  pulse 
code  modulation  systems  over  another. 

Implementation  Complexity  . The  complexity  of  the  transmitter, 
scanner,  and  the  receiver  all  depend  upon  the  scan  configuration.  The 
same  linear  array  of  N GaAs  laser  diodes  is  required  for  all  of  the  parallel 
scan  cases  and  the  serial  scan  case,  and  thus  will  be  equally  hard  to  imple- 
ment in  each  case,  whether  using  a single  substrate  or  N separate  diodes. 
The  mxn  array  required  for  the  scanning  array  probably  cannot  be  imple- 
mented on  a single  substrate,  and  would  be  several  orders  of  magnitude 
more  difficult  to  implement  with  individual  diodes.  As  mentioned  in  Chap- 
ter III,  the  optics  may  be  easier  to  implement  for  the  Overlapping  Beam 
Scan  configuration  than  for  the  other  scan  configurations.  For  the  remain- 
ing configurations,  the  optics  are  probably  of  the  same  degree  of  complexity 
from  case-to-case.  The  dimensions  of  the  rotating  prism  scanner  might 
be  changed,  and  the  rotation  rate  might  be  reduced  for  the  slow  parallel 
and  array  scan  cases,  but  otherwise  the  same  degree  of  complexity  is 
required  for  each  of  the  scan  configurations  considered. 

The  most  pronounced  differences  in  implementation  difficulty  from 
one  scan  configuration  to  another  occurs  in  the  receiver.  All  of  the 
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receivers  for  the  scan  configurations  above  require  N photodetectors, 
either  in  a linear  array,  or  in  an  mxn  array.  These  type  detectors  are 
available  and  pose  no  major  problems  other  than  the  non-uniformity  of 
responsivity  from  detector-ro-detector  which  occurs  in  both  the  linear  or 
rectangular  arrays.  Reviewing  the  analysis  results  of  Chapter  III  reveals 
that  the  Slow  Parallel  Scan  case  has  the  simplest  implementation;  it  requires 
no  delays  or  complicated  switching  network.  The  Series  Scan  case  was 
slightly  more  complicated,  requiring  one  delay  line  summer.  The  com- 
plexity increases  in  order  with  the  Array  Scan  and  Overlapping  Beams, 

Fast  Parallel  Scans,  and  Slow  Parallel  Scan-CP  Code;  requiring  m,  q,  and 
qN/2  delay  lines  summers,  respectively.  Hence  the  Slow  Parallel  Scan 
and  the  Serial  Scan  are  perhaps  the  easiest  scan  configurations  to  imple- 
ment when  using  pulse  code  temporal  modulation.  The  relaxation  of  the 
optics  complexities  for  the  Overlapping  Beams  configuration  may  overcome 
the  additional  required  receiver  complexities  (m  delay  lines)  for  the  scan 
system. 

Other  Considerations.  The  various  pulse  code  modulated  scan 
configurations  may  be  compared  in  several  other  respects,  such  as  the 
degree  of  raster  skew  from  the  perpendicular  to  the  aircraft  ground  track, 
and  the  number  of  averaged  transmitter- receiver  pair  responsivities.  In 
general,  the  faster  the  scan  rate,  the  smaller  the  skew  angle  will  be,  and 
visa-versa.  Therefore,  since  the  two  Fast  Parallel  Scanners  and  the 
Series  Scanner  sweep  at  the  same  rate  as  the  Single  Source  System,  they 
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will  have  identical  skew  angles.  The  Array  Scan  and  Overlapping  Beams 
configurations  are  scanned  at  1/m  times  the  rate  of  the  single  source,  and 
thus  have  larger  skew  angles.  The  Slow  Parallel  Scan  case  is  swept  at 
1/N  times  the  single  source  sweep  rate;  and  therefore,  has  the  largest  skew 
angle. 

As  was  pointed  out  in  Chapter  III,  the  Series  Scanner  and  the  two 
Fast  Parallel  Scanners  average  out  any  differences  in  the  N transmitter- 
receiver  responsivities,  thereby  ensuring  constant  reflectance  estimation 
performance  from  pixel-to-pixel.  The  Array  Scanner  and  Overlapping 
Beams  Scanner  only  average  n of  the  receiver  outputs,  thereby  allowing 
slightly  different  reflectance  estimation  performances  between  the  n lines 
of  each  scan.  The  Slow  Parallel  Scanner  does  not  average  the  outputs  of 
any  receivers,  and  thus  allows  for  even  greater  reflectance  estimation 
performance  from  line-to-line  of  the  scan.  It  should  be  pointed  out  that 
die  range  estimates  are  the  primary  data  required  for  target  cueing,  and 
the  reflectance  estimates  are  secondary  data.  Thus  tlie  loss  of  system 
performance  due  to  non- uniformity  of  reflectance  estimation  (performance 
depends  upon  the  degree  to  which  the  pattern  recognition  algorithms  require 
uniform  reflectance  estimation  (Performance,  and  may  not  be  too  severe  for 
the  differences  referred  to  here. 

It  follows  from  the  previous  discussion  that  the  Series  Scanner  has 
the  best  overall  PCM  coded  Scanner  Cliaracteristics.  It  has  the  same 
range  and  reflectance  (performance  as  tlie  other  configurations,  has  the 
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smallest  possible  raster  line  skew  angle,  has  uniform  reflectance  perform- 
ance from  line-to-line,  and  requires  only  one  delay  line  summer.  It  is 
also  the  simplest  method  of  implementing  the  linear  CP  code  modulation, 
which  has  the  most  desirable  autocorrelation  properties.  The  preceding 
comparisons  are  summarized  in  Table  I. 


Comparisons  of  the  Sinusoidal  Modulation  Class 

The  three  spatial- temporal  coding  combinations  in  the  sinusoid 
modulation  class,  "N  Sources  in  Parallel-Sinusoidal  Modulation /Multiple 
Sinusoidal  Modulation/OOK  Sinusoidal  Modulation",  will  now  be  compared 
with  each  other. 

Range  and  Reflectance  Performance.  The  range  performance  of  the 
three  systems  were  compared  in  Chapter  III  and  the  standard  deviation  (rms) 
of  the  range  error  of  the  systems  were  shown  to  be  related  by 


Vs  - ‘^r-M 

= ^ ^r-OOK 


(112) 


where  the  subscripts  S,  M,  and  OOK  denote  Single,  Multiple  and  On-Off- 
Keyed  Sinusoidal  Modulation,  respectively.  The  reflectance  performances 
were  also  compared,  yielding 


Op.s  = 1.1 

= /2  Op-oOK 


(113) 


where  <7p  is  the  reflectance  error  standard  deviation.  Thus  the  best 
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Comparison  of  PCM  Scan  Characteristics 


possible  range  and  reflectance  performance  is  achieved  by  the  Single  Sinu- 
soidal Modulation  System.  The  next  best  range  performance  is  (Atained 
by  the  OOK  Sinusoidal  Modulation  System.  The  reflectance  performance 
for  the  CXK  Sinusoid  case  is  only  29%  worse  than  that  of  the  Multiple  Sinu- 
soid case.  Since  the  reflectance  data  is  of  secondary  importance,  the  OOK 
Sinusc^dal  Modulation  system  is  the  second  best  sinusoidal  modulation  tech- 
nique for  estimating  range  and  performance. 

Maximum  Unambiguous  Range  . The  maximum  unambiguous  range 
for  the  three  sinusoidal  modulation  techniques  were  given  by 


^amb'S  ~ 2 fjjj 

Ramb-M  ' '^57 
Ramb-OOK  = 

2 

For  parameters  that  might  be  typical  of  laser  line  scan  systems,  assume 
f m = 100  MHz,  f = 1 MHz,  N = 10,  and  » 1 ^sec.  For  these  parameters, 
Eqs  (74),  (90),  and  (82)  yield  unambiguous  ranges  of  1.  5 meters,  150 
meters,  and  1500  meters  respectively.  Thus  the  later  two  modulation  tech- 
niques give  orders  of  magnitude  better  unambiguous  ranges  than  the  Single 
Sinusoid  case,  with  the  OOK  Sinusoidal  Modulation  case  yielding  the  best 
results. 


82 


l! 


Implementation  Complexity.  All  three  of  the  sinusoidal  modulation 
techniques  use  the  "N  Sources  in  Parallel"  scan  type  and  are  of  equal 
complexity.  The  modulation  of  the  drive  current  of  the  laser  diode  trans- 
mitters is  of  equal  complexity  for  the  Multiple  and  OOK  Sinusoidal  Modu- 
lation cases  and  is  only  marginally  more  complicated  than  that  of  the  Single 
Sinusoidal  Modulation  case.  The  receiver  configuration  complexity  increases 
in  order  with  the  Single,  Multiple,  and  OOK  Sinusoidal  Modulation  Techniques. 
However,  the  increased  complexity  of  the  OOK  Sinusoid  case  is  not  that 
severe  compared  to  the  Multiple  Sinusoid  case  and  is  certainly  worth  the 
large  increase  in  unambiguous  range  with  only  a /2  reduction  in  perform- 
ance compared  to  the  Single  Sinusoid  case.  Note  that  since  all  three  of 
these  sinusoidal  modulation  techniques  use  the  same  type  of  scan,  they  will 
have  identical  raster  line  skew  angles  and  reflectance  performance  behavior 
from  line  to  line. 

Comparisons  of  the  Non-Scanning  Class 

The  two  non-scanning  configurations,  "q-Starlng  Sources  and 
Detectors  - Pulse  Code  Modulation /Spatial  and  Temporal  Modulation"  may 
be  compared  by  inspection.  The  two  schemes  are  exactly  die  same  except 
the  latter  uses  spatial  coding  to  overcome  optics  deficiencies  in  order  to 
ensure  strict  temporal  ortiiogonality  between  adjacent  channels.  This 
reduces  the  effective  dwell  time,  and  tterefore,  the  efficiency  of  the  pulse 
code  modulation,  by  a factor  of  2/5.  Hence  the  standard  deviation  of  the 
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range  and  the  reflectance  errors  increase  by  a factor  of  1.  58  and  the  un- 
ambiguous range  is  decreased  by  a factor  of  1.  25  compared  to  the  staring 
scheme  which  does  not  employ  spatial  coding.  Thus  the  non-spatial- 
tempoial  coded  case  has  the  best  performance.  Note  that  the  use  of  the 
spatial  coding  merely  reduces  the  temporal  cross  talk  occurring  when  the 
illumination  from  one  channel  lies  in  the  field-of-view  of  another  channel. 
It  does  not  reduce  inter-pixel  spatial  cross  talk  due  to  the  illumination  of 
a given  channel  overspilling  the  pixel,  thereby  causing  "noisy"  range  and 
reflectance  estimates. 


Comparison  of  PCM  and  Sinusoidal  Modulations 

Range  and  Reflectance  Performance.  The  next  comparison  is  be- 
tween the  N Sources  in  Series,  Pulse  Code  Modulation  system  and  the 
N-Sources  in  Parallel,  Single  Sinusoidal  Modulation  System.  The  range 
performance  of  these  two  systems  can  be  compared  by  algebraically  com- 
bining Eqs  (80)  and  (86),  resulting  in 


®r-PCM  = T ^ ^r-S 

where  <^r-PCM  ^r-s  standard  deviations  of  the  range  error  of 

the  PCM  system  and  the  Single  Sinusoid  system,  respectively.  Likewise, 
the  reflectance  performance  can  be  compared  by  combining  Eqs  (81)  and 
(87),  yielding 


®P-PCM  0-87  a p.g 
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where  ®p-pcM  '^p-S  standard  deviations  of  the  reflectance 

error  for  the  PCM  systems  and  the  Single  Sinusoid  System,  respectively. 
Thus  the  reflectance  performance  of  the  PCM  system  is  always  better  than 
that  of  the  Single  Sinusoid  system.  From  Eq  (114),  it  is  easily  shown  that 
In  order  for  the  PCM  system  range  performance  to  be  greater  than  that  of 
the  Single  Sinusoid  system,  the  chip  duration  tQ  must  be  such  that 

‘o^^  = 0.64T„  016) 

where  is  the  period  of  the  modulation  sinusoid.  If  f^^  is  the  upper 
limit  of  the  frequency  response  of  the  laser  diode,  then  the  only  way  the 
PCM  system  can  be  made  to  outperform  the  single  sinusoid  system  is  to 
use  indirect  (electro-optic  or  spatial)  modulation.  If  fj^  is  the  upper 
limit  of  the  frequency  response  of  the  photodetector,  then  it  is  impossible 
for  tlie  PCM  system  to  outperform  the  Single  Sinusoid  System  in  range  esti- 
mation. 

Maximum  Unambiguous  Range.  The  maximum  unambiguous  range 
for  the  two  systems  are  given  by  Eqs  (83)  and  (74),  and  are 

cpto 

^'amb-PCM  ^ —7- 
^amb-S  “ 2ffn 

For  typical  system  parameters  of  pt^  = = 1 nsec  and  f(„  = 100  MHz, 

the  unambiguous  ranges  for  the  two  expressions  above  are  150  meters  and 
L 5 meters,  respectively.  Thus  it  is  seen  that  PCM  system  has  a superior 

I 

1 

I unambiguous  range  resolution  capability. 
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From  the  above  discussion,  it  can  be  concluded  that  the  Single 
Sinusoid  system  is  superior  if  the  short  range  ambiguity  is  acceptable. 

If  longer  unambiguous  ranges  are  required,  one  must  pick  an  acceptable 


range  ambiguity  and  then  use  Eq  (74)  to  calculate  the  required  f^.  The 
calculated  fj^,  and  the  smallest  chips  width  to  which  the  overall  system 
( is  capable  of  responding  to  must  then  be  used  in  Eqs  (80)  and  (86)  to 

determine  which  system  has  the  best  performance  for  the  required  unam- 
biguous range.  This  comparison  can  be  simplified  by  graphical  methods. 

Solving  Eq  (74)  for  f^^  and  substituting  into  Eq  (114)  yields 

" 4Ramb 

Solving  Eq  (117)  for  Ramb  setting  t^r-PCM  “ 'V-s  ^sults  in 
^annb  “ 4 *0 

which  may  be  plotted  as  shown  below  (Fig  28).  Now  if  the  intersection 
of  the  desired  unambiguous  range  and  available  state-of-the-art  t^  lies 
above  the  ttc  to/4  line,  the  PCM  system  will  give  the  best  performance; 
if  the  intersection  falls  below  the  line,  then  the  Single  Sinusoid  system  is 
the  best  choice. 

!l 

■J 


nc  / \ 

Ftg  28.  Graph  of  = ') 

Comparison  of  PCM  and  OOK  Sinusoidal  Modulations 

Range  and  Reflectance  Performance . The  next  comparison  is 
between  the  N Sources  in  Series,  Pulse  Code  Modulation  system  and  the 
N Sources  in  Parallel,  OOK  Sinusoidal  Modulation  system.  The  respec- 
tive range  performance  may  be  compared  by  combining  Eqs  (80)  and  (97), 
yielding  ^ 

®r-PCM  ^2/2  ^ Sr-OOK 

where  and  “^r-oOK  error  standard  deviations  for  tlie  PCM 

and  OOK  Sinusoidal  Modulated  systems,  respectively.  Likewise,  the 
reflectance  performances  may  be  compared  by  combining  Eqs  (81)  and 
(98),  giving 

VPCM  ‘^p-OOK 


i 


87 


w 


where  and  °p-OOK  reflectance  error  standard  deviations 

for  the  PCM  and  OOK  sinusoidal  modulated  systems,  respectively.  Thus 
the  PCM  system  will  always  have  the  better  reflectance  performance.  How 
ever,  the  PCM  system  will  have  the  better  range  performance  only  if 


0. 90 

= 0.9 


(121) 


As  stated  earlier,  if  the  frequency  response  of  the  system  is  limited  by 
the  laser  diode,  then  indirect  modulation  must  be  used  in  order  for  the  PCM 
system  to  have  better  range  performance  than  the  OOK  Sinusoid  system. 

If  the  system  frequency  response  is  limited  by  the  detector,  then  the  OCX 
Sinusoid  system  will  always  have  the  better  performance. 

Maximum  Unambiguous  Range  . The  maximum  unambiguous  range 
of  the  N Sources  in  Series,  PCM  Modulation  system  and  that  of  the  N 
Sources  in  Parallel,  OOK  Sinusoidal  Modulation  system  may  be  compared 
by  combining  Eqs  (82)  and  (83),  yielding 


^mb-OOK  ^ ^amb-PCM 


(122) 


Note  that  the  restriction  on  t^  given  by  Eq  (102)  must  be  met  in  order  for 
Eq  (122)  to  be  valid.  So  long  as  that  restriction  is  satisfied,  there  is  no 
advantage  to  using  the  PCM  system  to  resolve  range.  Thus,  based  on 
the  above  analysis  of  all  the  scanning  systems  using  N sources,  the  OOK 
Sinusoidal  Modulation  system  has  the  best  range  performance  if  large 
range  ambiguity  is  required.  If  large  range  ambiguities  are  not  required. 
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tfien  the  Single  Sinusoidal  Modulation  system  has  the  best  range  perform- 
ance. 

Implementation  Complexity  . Both  the  Series  Scan,  PCM  system 
and  the  Parallel  Scan,  OOK  Sinusoidal  Modulation  system  use  the  same 
linear  array  of  laser  diodes,  the  same  optics  and  except  for  perhaps  minor 
changes  in  dimensions  or  rotational  rate,  the  same  scanner.  The  receiver- 
processor  for  the  OOK  Sinusoid  system  is  more  complicated  than  that  of 
the  PCM  receiver  due  to  the  number  of  PLL,  mixers,  and  filters  required. 
However,  it  is  much  more  practical  to  realize  than  some  of  the  Fast  Para- 
llel Scan  system  described  earlier,  and  is  certainly  not  unreasonable. 
Sinusoids  are  also  easier  to  generate  than  nanosecond  switching  networks 
and  delays. 

Other  Considerations.  The  PCM  system  above  employs  the  serial 
scanner  and  thus  has  uniform  reflectance  estimation  performance  from  line- 
to-line  as  well  as  a small  raster  line  skew  angle.  The  OOK  Sinusoidal 
Modulation  system,  on  the  other  hand,  employs  the  parallel  scanner.  Thus 
the  OOK  Sinusoidal  Modulation  system  does  not  average  any  of  the  receiver 
output;  and  therefore,  may  not  necessarily  exhibit  uniform  reflectance 
estimation  performance  from  line-to-line  of  the  raster  scan.  However, 
depending  upon  the  requirements  of  the  pattern  recognition  algorithm, 
this  may  or  may  not  be  a serious  problem.  The  OOK  Sinusoidal  Modu- 
lation system  also  has  a larger  raster  line  skew  angle  than  the  PCM  system. 
However,  this  is  not  a major  concern  since  the  raster  lines  may  be  de- 
skewed by  die  pattern  recognition  processor. 
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It  is  important  that  one  other  point  be  recognized  here.  The  above 
analysis  comparing  the  PCM  Serial  Scan  system  to  the  OOK  Sinusoidal 
Modulation  system  could  be  generalized  to  include  the  other  scan  configu- 
rations that  employ  PCM  with  few  changes.  However,  if  the  Overlapping 
Beams  scan  configuration  is  not  restricted  to  N sources,  then  the  general 
comparisons  above  will  not  apply.  Through  the  brute  force  technique  of 
adding  more  sources;  and  consequently,  more  energy  per  pixel;  the  fre- 
quency response  deficiency  of  the  PCM  coding  scheme  can  be  overcome. 
Thus,  the  Overlapping  Beams,  PCM  coding  scheme  can  be  made  to  out- 
perform the  Single  or  die  OOK  Sinusoidal  Modulation  coding  schemes. 

This  system  also  has  the  advantages  of  averaging  N receiver  outputs  per 
scan  line,  a reduced  raster  line  skew  angle  and  possibly  simpler  trans- 
mitter optics. 

Comparison  of  q-PCM  Channels  and  N-OOK  Sinusoidal  Modulation  Channels 

Range  and  Reflectance  Performance.  The  final  two  systems  to  be 
compared  are  the  N Sources  in  Parallel,  OOK  Sinusoidal  Modulation  sys- 
tem and  the  q-Staring  Sources  and  Detectors,  PCM  system.  The  range 
performance  of  the  two  systems  may  be  compared  by  combining  Eqs  (97) 
and  (103),  resulting  in 

Sr-N-OOK  = • Vq-PCM  (123) 
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where  <^i-n-ook  '^r-q-PCM  range  estimation  errors 

for  the  N"OCK  source  system  and  tlie  q-PCM  source  system,  respectively. 
The  reflectance  performance  of  the  two  systems  may  be  compared  by  com- 
bining Eqs  (98)  and  (104),  yielding 

Vn-OOK  = * ^p-q-PCM  (^24) 


where  ^p-jsi-ocK  '^p-q-PCM  reflectance  estimation 

errors  for  the  N-OOK  source  system  and  the  q-PCM  source  system,  res- 
pectively. The  above  results  are  hard  to  interpret  because  of  the  many 
independent  variables.  However,  the  range  estimation  performance  of 
tlie  N-source  system  can  only  exceed  that  of  the  q-source  system  if 


0.9 


7 


m^o 


< 1 


(125) 


Solving  Eq  (125)  for  N gives 


026) 


Two  examples  will  give  a better  insight  as  to  the  relationship  between 

®r-N-OOK  Vq-PCM-  ^ ^o  " ^/^m’  N >33 

in  order  for  the  N-OOK  source  system  to  have  tlie  smaller  rms  range 

error.  If  q = 10,000  and  tp  = 10/fj^ , then  N need  only  be  81  to  out- 
perform the  10,000  element  staring  system.  This  indicates  that  the  im- 
provement in  the  range  performance,  which  is  proportional  to  the  square 
root  of  the  number  of  sources,  is  not  worth  the  additional  complexity 
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Involved  for  the  q-source  staring  system.  These  results  are  general 
enough  tliat  they  can  be  expanded  to  the  Overlapping  Beams,  PCM  case 
also,  wliere  the  number  of  overlapping  sources  is  allowed  to  be  q,  q > N. 
Thus,  if  a laser  line  scanner  is  required  to  have  a long  unambiguous  range, 
then  the  N Sources  in  Parallel,  OOK  Sinusoidal  Modulation  combination  of 
spatial  and  temporal  coding  will  yield  the  best  overall  range  and  reflectance 
estimation  performance.  The  comparisons  results  of  this  chapter  are  sum- 
marized in  Table  11. 

Table  11 

Spatial -Temporal  Systems  Cliaracteristics 

Range 


System 

^r 

Ambiguity 

P Interval 

Optics 

Complexity 

Serial  Scan,  PCM 

TT  - . 

•y  *^nito'^r-s 

Vpcm  ^\i/2 

Overlapping  Beams,  PCM 

TWo'^r-S 

°p-PCM  c'^d/2 

Least 

ftirallel  Scan,  Single 
Sinusoidal  Modulation 

‘^r-S 

I.15ap_ppj^  c/2f,T^ 

Severe 

l^rallel  Scan,  OOK 
Sinusoidal  Modulation 

l.bSOp-pcM  cT^/2 

As  an  example  of  the  range  performance  which  might  be  expected 
with  a multiple  source  GaAs  laser  line  scan  system,  assume  tlie  following 
system  parameters:  N = 10,  = 1 watt.  Eg  q 35^^  = 0.  07  watts/cm^-p  m, 

p=0. 1,  6=0°,  tj.  = tj  « I,  Ai.  = 100cm*,  Aj  = 900  cm’’ (~lft  x 1ft),  r = 

2.  5 x 10*  cm  (~820  ft),  A \ = 0.01  pm,  f^i  Mllz,  and  Tj  = 1 psec. 
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Using  Eq  (77),  (78)  and  (97),  recognizing  that  = ^Pj^,  the  rms  range 
error  for  the  OOK  Sinusoidal  Modulation  system  is  '^r-ocK  “ 2.  3 x 10“\:m, 
certainly  good  enough  target  for  cueing.  Note  that  this  is  the  best  per- 
formance possible,  theoretically.  Actual  attainable  rms  error  may  be 
larger.  Also,  this  example  was  calculated  assuming  no  losses  due  to  the 
atmosphere  or  optics,  and  neglected  noise  due  to  atmospheric  backscatter. 
For  operation  in  fog  or  haze  the  performance  would  degrade  appreciably. 
Thus  die  use  of  several  sources  for  clear  weather  operation  may  seem  like 
overkill,  but  may  be  quite  necessary  of  inclement  weather  operation.  Also, 
note  that  the  rms  range  error  varies  as  fj^  and  as  (NP^j)"^.  Thus  it  is 
more  important  to  use  sources  with  good  high  frequency  response  than  with 
high  output  powers,  and  the  number  of  sources  N required  will  depend  upon 
fmt  Pjj»  ®nd  the  performance  required. 


Comparisons  Results:  rms  Power  Limited  Sources 

All  of  the  preceding  analysis  assumed  the  GaAs  laser  diodes  were 
peak  power  limited.  If  this  assumption  were  removed,  the  above  results 
would  change  only  slightly.  Since  Eq  (39)  indicates  that  the  reflectance 
performance  depends  only  upon  the  post  detection  signal  energy,  removal 
of  the  peak  power  limited  restriction  means  that  all  of  the  spatial-temporal 
coding  combinations  will  have  identical  reflectance  performances.  As  for 
the  range  performance,  removal  of  die  restriction  enables  the  performance 
of  the  q-Staring  Sources  and  Detectors,  Spatial  and  Temporal  Modulation 


W 


. . V 


I 


case  to  equal  that  of  the  q-Staring  Sources  and  Ctetectors,  PCM  case.  It 
also  allows  die  reflectance  performance  of  all  of  the  sinusoidal  cases  to 
be  equal.  However,  it  can  be  shown  tliat  if  Eqs  (80)  and  (86)  are  normal* 
ized  by  Eg.p^f^  and  , respectively,  and  the  algebraically 

combined,  the  result  is 


1 1 

‘^r-PCM  “ 72  ^m*^o  S:-s 


Therefore,  it  is  required  that 


/2 
^ tt 


027) 


028) 


in  order  for  the  PCM  system  to  outperform  the  sinusoidal  system.  There- 
fore, even  if  the  peak  power  limited  restriction  is  removed,  the  best  range 
performance  is  still  achieved  by  using  sinusoidal  modulation  or  OOK  sinu- 
soidal modulation. 
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V.  Conclusions  and  Recommendations 


This  Chapter  contains  the  general  conclusions  based  upon  the  anal- 
ysis and  comparisons  of  the  previous  chapters,  and  lists  several  recom- 
mendations for  further  investigation. 

Conclusions 

The  basic  purpose  of  this  study  was  to  find  the  combination  of  scan 
configuration,  for  N GaAs  laser  diode  sources,  and  intensity  modulation 
temporal  coding  which  would  produce  optimum  range  and  reflectance  esti- 
mation performance.  Toward  this  end,  the  statistical  additive  noise  model 
for  a direct  (optical  intensity)  detector  was  developed.  Estimation  theory, 
based  upon  the  conditional  Poisson  statistics  of  the  detector,  was  used  to 
develop  the  joint  maximum-likelihood  (ML)  estimator  (mathematical  model) 
for  both  range  and  reflectance.  The  realization  of  the  joint  estimator  was  a 
correlator  which  required  the  reflectance  estimate  times  a replica  of  the 
transmitted  signal  (delayed  by  the  time  estimate)  for  an  input.  The  joint 
estimator  could  not  be  easily  realized.  However,  it  was  shown  that  a real- 
ization of  the  range  estimator  is  a matched  filter,  and  a suboptimum  reali- 
zation of  the  reflectance  estimator  is  the  measurement  of  the  detected  sig- 
nal energy,  given  the  range  estimate.  The  Cramer- Rao  lower  bound  was 
calculated  for  both  range  and  reflectance  estimates.  The  variance  of  the 
range  error  is  inversely  proportional  to  the  mean  square  bandwidth  of  the 
modulating  signal,  and  the  variances  of  both  the  range  and  the  reflectance 


errors  are  inversely  proportional  to  the  detected  signal  energy.  Various 
combinations  of  spatial  and  temporal  coding  of  the  N laser  sources  to 
separate  the  return  signal  energy  from  each  ground  resolution  cell  while 
increasing  the  returned  energy  per  cell  were  evalua«'ed.  This  was  done  by 
analyzings  large  variety  of  spatial-temporal  and  hybrid  spatial- temporal 
coding  schemes  and  then  comparing  the  performance  results.  Unfortunate- 
ly, one  technique  is  not  always  superior;  the  best  scanning-coding  combi- 
nation to  use  depends  upon  the  requirements  of  the  pattern  recognition- 
algorithm.  Nonetheless,  several  general  conclusions  maybe  made  and 
will  be  discussed  below. 

The  suboptimum  matched  filter/peak  detector  realization  of  the 
joint  estimator  is  suboptimum  because  the  round  trip  delay  estimate  is 
required  in  order  to  obtain  the  reflectance  estimate.  Therefore,  the  re- 
flectance measurement  is  dependent  upon  die  delay  measurement.  However, 
this  appears  to  be  the  best  possible  processor  configuration  which  can  be 
constructed  using  current  technology.  The  performance  predicted  by  the 
Cramer  Rao  lower  bound  on  the  mean  square  estimation  error  is  the  smal- 
lest error  limit  theoretically  possible.  Regardless  of  the  processor  used, 
a performance  better  than  the  Cramer  Rao  lower  bound  can  not  be  achieved. 
In  practice,  this  lower  bound  may  not  be  reached  easily.  Yet  the  Cramer 
Rao  lower  performance  bound  serves  as  a standard,  or  bench  mark 
against  which  the  performance  of  real  systems  may  be  compared.  And 
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since  the  performance  lower  bounds  calculated  for  the  systems  suggested 
In  this  thesis  are  well  within  an  acceptable  range,  it  is  practical  to  try 
to  build  a laser  line  scan  sensor  using  N GaAs  semiconductor  laser 
sources. 

If  the  main  goal  of  die  laser  line  scan  system  design  is  to  achieve 
the  smallest  possible  range  and  reflectance  errors,  and  if  the  short 
ambiguity  interval  and  lack  of  uniform  reflectance  estimation  perform- 
ance from  Une-to-line  is  not  a problem,  then  the  parallel  scan  configuration 
(Fig  12,  page  43)  using  sinusoidal  intensity  modulation  is  the  best  choice. 
However,  if  longer  ambiguity  intervals  are  required,  then  the  same  para- 
llel scan  configuration  using  On-Off-Keyed  sinusoidal  modulation  appears 
to  be  the  preferred  technique.  If  in  addition  to  long  ambiguity  intervals, 
uniform  reflectance  estimation  performance  from  line-to-line  is  required, 
then  either  the  pulse  code  modulated  serial  scan  configuration  (Fig  15, 
page  50)  or  the  pulse  code  modulated  overlapping  beams  scan  configuration 
I (Fig  26,  page  72)  may  be  used.  In  the  latter  case,  the  overlapping  beams 

configuration  may  have  slight  reflectance  performance  variatiOTs  from 
line-to-line,  but  may  employ  simpler  transmitter  optics.  All  of  the  other 
techniques  require  that  the  N individual  transmitter  fields-of-view  do  not 
overlap,  i.  e. , are  spatially  orthogonal,  and  all  of  the  above  techniques 
require  the  receiver  fields-of-view  to  be  spatially  orthogonal.  It  is  con- 
eluded  that  the  most  cost  effective  system  with  the  best  overall  range, 
reflectance,  and  maximum  unambiguous  range  characteristics  available 
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using  existing  technology  is  the  N Sources  in  Parallel,  OOK  Sinusoidal 


Modulation  scheme.  The  sinusoidal  modulation  frequency  f should  be 

m 

as  high  as  possible,  and  will  be  limited  by  the  upper  frequency  response 
of  the  source  or  the  photodetector,  whichever  is  lower.  Narrow  (fast 
lisetime)  optical  pulses  are  not  required,  thus  wide  bandwidths  are  not 
required.  The  detector  need  only  be  fast;  it  does  not  need  to  have  DC  to 
RF  frequency  response.  High  peak  power  GaAs  laser  diodes  are  not  essen- 
tial, but  high  peak  powers  combined  with  high  frequency  responses  ensure 

better  performance  in  inclement  weather  and  reduces  the  number  of  sources 

\ 

required. 

Recommendations 

Based  on  the  conclusions  above,  the  first  recommendation  for 
further  work  in  this  area  is  for  a detailed  study  to  examine  existing  or 
proposed  cueing  algorithms  in  order  to  determine  what  ambiguous  range 
interval  and  reflectance  performance  criteria  are  required.  This  would 
then  be  used  to  determine  which  of  the  four  techniques  above  should  be 
emphasized  die  most  in  future  research. 

Another  area  requiring  further  investigation  is  that  of  spatial  fil- 
tering. The  scanning  process  was  modeled  as  if  the  illuminating  beam 
moves  discretely  from  pixel  to  pixel  with  no  smearing  due  to  the  beam 
movements  or  forward  motion  of  the  aircraft.  There  are  some  indica- 
tions that  the  spatial  filtering  may  be  significant:  therefore,  the  area 
warrants  additional  investigation. 
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In  order  to  limit  the  scope  of  the  study,  the  only  pulse  modulation 
codes  considered  were  Barker  codes,  pseudo-noise  (PN)  codes,  and 
Totally  Orthogonal  Complimentary  Pair  (CP)  codes.  There  are  other 
coding  schemes  with  long  codes  and  bi-level  autocorrelation.  This  is  an 
area  which  deserves  more  attention. 

This  analysis  assumed  that  tlie  Fields-of-View  (FOV)  of  the  trans- 
mitter optics  and  the  receiver  optics  could  be  matched,  and  could  be  made 
spatially  orthogonal  to  adjacent  transmitter-receiver  pair  FOV.  Since  the 
GaAs  laser  diode  output  beam  has  a rather  large  divergence,  this  may  not 
be  true.  Therefore,  the  feasibility  of  implementing  the  optics  for  each  of 
the  recommended  scan  techniques  above  should  be  evaluated.  Also,  the 
problem  of  summing  the  outputs  of  N laser  diodes  into  a single  narrow 
beam  deserves  further  investigation.  Further  study  is  needed  to  deter- 
mine the  central  beam  characteristics  of  GaAs  source  and  how  to  control 
the  beam  characteristics  using  both  diode  fabrication  and  optics. 

The  final  recommendation  is  for  further  work  in  component  tech- 
nology. A concentrated  effort  to  develop  linear  arrays  of  laser  diodes, 
and  of  photodetectors,  with  uniform  performance  from  element  to  element 
is  desirable.  An  initial  goal  of  five  or  ten  GaAs  laser  diodes  in  a linear 
array  might  be  reasonable. 
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Appendix  A | 

i 

Additive  Noise  Model  for  Direct  Detection 

The  purpose  of  this  appendix  is  to  develop  the  additive  noise  model 
for  a direct  detector  shown  in  Chapter  II  (Fig  1).  The  total  detector  out- 
put current  is 

The  mean  of  the  filter  output  will  be  the  expected  value  of  i^(t)  convolved  . 
with  the  impulse  response  of  the  filter.  Therefore, 

E[y(t)]  = E[i^(t)]  * h(t)  030) 

In  order  to  evaluate  E [ijj(t)] , one  may  use  the  fact  that 

l(t)  = q iimi.  031)  1 

dt 

where  N is  the  total  number  of  counts  (detected  photons)  over  the  obser- 
vation interval.  Using  Eq  031).  one  finds  that  in  general  (Ref  4 : 258-260) 

ECi(t)]  = qX(t)  032) 

and 

Ri(t.f)  = q»X(t)6(t-t’)  + q»X(t)X(t’)  033)  | 

Eq  (132)  and  Eq  033)  imply  that  the  combination  of  dark  current  plus 
background  noise  in(t)  = ijj(t)  + iQ(t)  may  be  modeled  as 


ECin(t)]  = q\i 

Ri^[t,t']  = 6(t-t’) 


(134) 

(135) 


where  ' 1^  was  mentioned  in  Chapter  II  that  the  thermal 

noise  has  first  and  second  moments 

ECLj^(t)]  = 0 (136) 

and 

KT 

= 2 — 6(t-t’)  (137) 

Using  Eqs  (130),  (132),  (133),  (134),  and  (135),  becomes 

E[id(t)]  = qXg(t)  + qXj^  038) 

Writing  Eq  (130)  in  integral  form  and  substituting  Eq  (138)  yields 


E Cy(t)]  = IZ  h(t-a)  q (a)  + da  039) 

The  autocorrelation  of  the  filter  output  will  be  the  autocorrelation 
of  ijj(t)  ♦ h(t)  and  may  be  expressed  as 

Ry(^  t’)  = n da  dp  h(t-a)  h(t’-p)  E Cid(a)i^(P)]  (140) 

Using  E[X]  = E[X]ECY]  for  X and  Y statistically  independent,  and 

substituting  Eq  (129),  and  Eqs  (131)  thru  (136)  and  rearranging  gives 

KT 

E Cid(a)id(P)]  = {q**Cya)+Xn]  + 2-^  } 6(a  - p) 

+ q**  lXg(a)Xg(p)  + X^[Xg(a)  + X Jp)]  + X»^  ) (141) 

The  covariance  of  the  filter  output  is  given  by 
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1 

Cy(t,f)  = Ry(t.t’)  - E [y(t)]  E [y(t’)]  (142) 

Substituting  Eqs  (138)  thru  (141)  into  Eq  (142)  and  simplifying  yields 
Cy(t,n  = Hda  Jl*d3  {qnXg(a)+X^]+2^  3 6(a-3)  043) 

Therefore,  zero  mean  additive  noise  model  may  be  represented  by  letting 
the  input  be  E [i(j(t)]  and  the  additive  noise  be  n(t)  with  Rj^(t,  t')  = C^(t,  t’). 

The  schematic  of  the  resulting  noise  model  is  shown  in  Chapter  n (Fig  1). 
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